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ABSTRACT 

We report observations of Faraday rotation measures for a sample of 191 extragalactic radio jets 
observed within the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 
program. Multifrequency Very Long Baseline Array (VLBA) observations were carried out over twelve 
epochs in 2006 at four frequencies between 8 and 15 GHz. We detect parsec-scale Faraday rotation 
measures in 149 sources and find the quasars to have larger rotation measures on average than BL Lac 
objects. The median core rotation measures are significantly higher than in the jet components. This 
is especially true for quasars where we detect a significant negative correlation between the magnitude 
of the rotation measure and the de-projected distance from the core. We perform detailed simulations 
of the observational errors of total intensity, polarization and Faraday rotation, and concentrate on the 
errors of transverse Faraday rotation measure gradients in unresolved jets. Our simulations show that 
the finite image restoring beam size has a significant effect on the observed rotation measure gradients, 
and spurious gradients can occur due to noise in the data if the jet is less than two beams wide in 
polarization. We detect significant transverse rotation measure gradients in four sources (0923-1-392, 
1226+023, 2230-1-114 and 2251-1-158). In 1226-1-023 the rotation measure is for the first time seen to 
change sign from positive to negative over the transverse cuts, which supports the presence of a helical 
magnetic field in the jet. In this source we also detect variations in the jet rotation measure over a 
time scale of three months, which are difficult to explain with external Faraday screens and suggest 
internal Faraday rotation. By comparing fractional polarization changes in jet components between 
the four frequency bands to depolarization models we find that an external purely random Faraday 
screen viewed through only a few lines of sight can explain most of our polarization observations but 
in some sources, such as 1226-1-023 and 2251-1-158, internal Faraday rotation is needed. 
Subject headings: BL Lacertae objects: general - galaxies: active - galaxies: jets - polarization — 
quasars: general - radio continuum: galaxies 



1. INTRODUCTION 
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Polarimetric observations of active galactic nuclei 
(AGN) jets enable studies of the magnetic field struc- 
ture in the outflows. If the jets are launched from the 
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rotating black hole or accretion disk, it is natural to ex- 
pect that the magnetic field structure in the jets is he- 
lical (e.g., 'Blandford fc Znaiek'l l977l [Mei er et al."2001^ 
Vlahakis fc K5nigl 2004; McKi nney fc N arayan 2007). 
Un the other hand, it is not known whether the helical 
structure persists parsecs down from the central engine 
or if it becomes tangled due to re-collin iation shocks or 
intera ction with external medium (e.g., iMarscher et al. 



2008|. Very Long Baseline Interferometry (VLBij can 
be used to study the electric vector orientation in the 
parsec-scale jets of AGN. In the optically thin part of 
the jet, the magnetic field orientation is perpendicular 
to the electric vector position angles (EVPAs). Thus, 
observations of EVPAs parallel to the jet direction have 
result ed in claims of toroida lly dominated magnetic fields 
(e.g., |Gabuzda et al.||2004|). One should note that rela- 



tivistic ettects make the situation more complicated and 
when viewed at small angles a toroidally dominated mag- 
netic field can appea r poloidal in the observer's frame 
( Lyutikov et al.|[2005 ). Alternatively the observed mag- 
netic field orientation can be accounted for by shocks 
compressing the magnetic field perpendicular to the jet 
( [Laing 1980] [Hughes et aI]|T989j ). 

Polarized waves are atiected by Faraday rotation when 
propagating through non-r elativistic p lasma within or 
external to the source (e.g., [Rum|p66| ). This effect can 
both diminish the observed degree ol polarization and 
rotate the intrinsic EVPAs so that in order to study the 
intrinsic magnetic field orientation in the jets, the effect 
must be removed. In the case of external rotation, the 
effect can be described by a linear dependence between 
the observed EVPA (xobs) and wavelength squared (A^) 

by 



Xobs = XO + 



y"neB-dl = xo + RMA2, (1) 



where xo is the intrinsic EVPA and RM is the rotation 
measure, related to the electron density rie and the mag- 
netic field component B parallel to the line of sight. The 
constant in the equation consists of the charge of the 
electron e, vacuum permittivity eg, mass of the electron 
m, and speed of light c. The RM can thus be estimated 
by observing the EVPA at several frequencies. 

If the Faraday rotation is internal to the jet it means 
that either the thermal plasma causing the rotation is 
intermixed with the emitting plasma, or the relativis- 
tic particle spectrum extends to low energies. For total 
rotations larger than 45° internal Fara day rotation is ex- 
pected to cause s evere depolarization (|Burn|1966 ) which 
is not often seen f Zaval a fc Taylor||200l ), but tor smaller 
total rotations the internal rotation can appear linear and 
follow Eq. [l] External Faraday rotation could be caused 
by a screen very close to the jet itself where it can also 
interact with the jet (e.g., a sheath) or by a more dis- 
tant screen such as the broad or narrow line regions, or 
even intergalactic and Galactic plasma. Distinguishing 
between the different alternatives can be very difficult, 
especially in the case of small r otations if addition al in- 
formation is not available (e.g., Homan et al.||2009 ). 

Over the past few decades there have been numerous 
studies of Faraday rotation in parsec-scale jets a ssoci- 
ate d with active galaxi es. One of the largest is by|Tay- 
lor| ( |l998l [2000) and ,Zavala fc Taylor] ( [20^ [2003| [200% 



who report Faraday rotation measures (RMs) for a sam- 
ple of 40 AGN. They find the typical absolute core RMs 
to be in a range of 500 to several thousand rad in 
the observer's frame. Additionally they report variabil- 
ity in several RMs over a time span of months to years, 
ruling out the narrow line region as the origin for the 
Faraday rotation. Instead they suggest that the Faraday 
rotation is caused by a screen close to the jet. Similar 
conclu sions about the screen were drawn by |Asada et al.| 
( 2002 1 who detected a transverse rotation measure gra- 
dient in the jet of 3C 273. They interpreted the gradient 
as a signature of a helical magnetic field in the sheath 
surrounding the jet. Several other claims of transverse 
gradients hav e been published (e.g. , Gabuzda et al. 12004 



Asada et aLl|2008b[ [Gomez et aLl|2008| [Mahmud et aL 
2009| |Croke et al' 2010p but due to some uncertainties 

regarding transversely unresolved jets the issue remains 
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controversial ( jTaylor fc Zavala||2010[ ). 

Due to the complex nature of the sources and their 
surroundings, the situation is often more complicated. In 
addition to a helical field, interactions with the surround- 
ing intergalactic medium will cause distinct RM struc- 
tures ( Gomez et al. 2008) . Additionally, beam effects 
can severely complicate tne interpretation of the map s 
as shown by simulations ( [Broderick fc McKinney|[2010 ). 
Starting from g eneral relativistic magnetohyd rodynami- 
cal simulations, Broderick fc McKinney (2010) created a 
canonical jet model and calculated RM maps which they 
convolved with different image restoring beam sizes to 
create unresolved and resolved jets. They showed that 
within one beam width from the optically thick core, any 
gradient seen in the RM map is generally unreliable, and 
only at resolution obtained with 43 GHz VLBA obser- 
vations are in agreement with expected values, although 
still suppressed in magnitude. In the case of an opti- 
cally thin jet, it could be possible to detect a gradient 
if the jet is surrounded by a helical field even if the jet 
is unresolved (above 8 GHz VLBA resolution), but the 
magnitude of the gradient may be suppressed due to the 
beam effects. 

In this paper we study the statistical properties of 
Faraday rotation in AGN by using a large sample of ob- 
jects which are part of the MOJAVE (Monit oring of Jets 
in AGN with VLBA Experiments) survey (Lister et al. 
[2009a). Our goal was to create a set of RM maps in which 
all potential sources of error in the data processing have 
been accounted for. Therefore we have performed exten- 
sive simulations of the errors in polarization and Faraday 
rotation maps and have assessed when a RM gradient can 
be called significant. These simulations show that the fi- 
nite beam size of VLBI observations has a large effect on 
the observed Faraday rotation, and caution needs to be 
taken when interpreting the maps. Neighboring pixels in 
beam-convolved images are not independent for a typical 
VLBI pixel size and the RM maps generally consist only 
of a few independent measurements. 

We describe our observations and the detailed data 
analysis process in Section [2l The results of our statis- 
tical study are reported in section [3| We discuss our 
results in light of depolarization models, observed RM 
gradients and time variability in Section [4j Our conclu- 
sions are summarized in Section[5] In Appendix A we dis- 
cuss the effects of relative image alignment on RM maps, 
and in Appendices B and C we discuss observational er- 
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rors in polarization and RM images. Throughout the 
paper we use a cosmology whe re Hp = 71 kms~^Mp c~^, 
= 0.3, and f^A = 0.7 (e.g., |Komatsu et al.||2009 |. 



2. OBSERVATIONS AND DATA REDUCTION 

Our sample consists of 1 91 AGN observed w ithin the 
MO J AVE VLB A survey ( [Lister et aT]|2009a[ ). It in- 
cludes 134 sources of the complete ttux density- limited 
MOJAVE-1 sample, for which we monitor total intensity 
and polarization changes of 135 AGN jets above decli- 
nation —20° which have exceeded 15 GHz flux density 
1.5 Jy (2Jy at 5 < 0°) at any epoch between 1994 and 
2002. The rest of the sources belong to the MOJAVE- 
2 sampl^ which includ es sources from the 2 cm survey 
( [Kellermann et al.||2004| , gamma-ray blazars, and other 
sources with unusual jet properties. 

The sources were observed with the VLBA in 2006 
over 12 epochs with about monthly separation, each 
epoch containing 18 sources (except for epoch 2006-Feb- 
12 which included only 14 sources and epoch 2006-Apr-28 
which included 17 sources). The observations were made 
in dual polarization mode using frequencies centered at 
8.104, 8.424 (X-band), 12.119 and 15.369 GHz (U-band). 
This setup was chosen because VLBA observes the gain 
centered at 8.4 GHz, and 8.1GHz was chosen as the low 
frequency end of the X-band. The bandwidths were 16 
and 32 MHz for the X and U-bands, respectively. The 
observations were recorded with a bit rate of 128 Mbits 
s^^. In the X-bands the observations consist of 2 IFs 
in both frequencies and in the U-bands 4 IFs. All ten 
VLBA antennas were observing except at epoch 2006- 
Aug-09 when Pie Town was not included. The sources 
and their observing epochs are listed in Table [T] A total 
of twenty sources were observed twice during the year. 

2.1. Data reduction 

The initial data reduction and calibration were per- 
formed following the standard procedures described in 
the AIPS cookboolr] All the frequency bands were 
treated separately throughout the data reduction pro- 
cess. The imaging and self-calibration were don e in a 
largely aut omated way using the D ifmap packag e ( Shep 



herd 
loFT 



T997| . For more details see 



le standard d ata reduction and imaging process and 
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Lister & Homan (2005) for the calibration of the polar- 
ization data. 

All the maps were modelfit with circular or elliptical 
Gaussian components using the standard procedure in 
the Difmap package. The 15 GHz maps were modelfit 
already as a part of the MOJAVE survey (Lister et al. 
12009b). Since one of our goals was to use the optically 
thin components in the jets to align our images, we used 
these 15 GHz modelfits as a starting point for the other 
bands and modified the fit if needed. 

As the {u,v) plane coverage differs in the bands with 
higher frequency maps resolving smaller structures, we 
can get spurious features in the rotation measure maps, 
especially near the core region where there can be many 
components blending within the beam at lower frequen- 
cies. Therefore, in order to have comparable {u,v) cover- 
age in all the bands, we flagged the long baselines from 



^ http:/ /www. physics. purdue.edu/astro/MOJAVE/allsourccs.shtml 
^ http://www.aips.nrao.edu 



the 15 and 12 GHz maps and short baselines from the 
8 GHz maps. The resulting typical {u,v) range in our 
data is 7.3 - 231 mega-A. We tested on individual sources 
the difference between flagging the baselines compared to 
tapering and found the differences to be minimal; there 
was no difference in the final RM map area, and the dif- 
ferences in RM values were a fraction of the error bars. 
Additionally, we restored all the maps to the beam size 
of our lowest frequency (8.1GHz). All these steps were 
carried out in Difmap after the initial data reduction and 
self-calibration which were done using the full {u,v) data. 

2.2. Absolute EVPA calibration 

The absolute electric vector position angle (EVPA) off- 
set is an instrumental quantity that must be determined 
and applied to every VLBA polarization observation. To 
calibrate the EVPAs of our data, we used Very Large 
Array (VLA), University of Michigan Radio Astronomy 
Observatory (UMRAO) and 15 GHz VLBA data, and in- 
strumental leakage term (D-term) phases. The 15 GHz 
observations were previously calibrated as part of the 
MOJAVE project using the D-term c alibration method 
( [Gomez et al.|2002[pster fc Homan]20 05). Therefore we 
only had to calibrate the 8 and 12 GHz bands. For five 
epochs we were able to use the VLA/ VLB A polarization 
calibration database]^ to find polarization observations 
within a week of our epoch and including one or two of 
our sources. For those epochs, we also calculated the 
distribution of differences between the calibrated 15 GHz 
EVPAs and other bands, and UMRAO 8 and 15 GHz 
EVPAs versus our 8 and 15 GHz EVPAs. Usually these 
difference histograms showed a peak at an angle similar 
to that determined from the VLA observations. The typ- 
ical errors in the VLA EVPAs range from 1 to 3 degrees 
and in the UMRAO data from 1 to 10 degrees but these 
cancel when multiple sources are used. 

By using these five epochs, we were able to find D-term 
phases on various antennas that were stable enough over 
the 12 month period to enable the use of D-term phases 
in the calibration of the EVPAs of the remaining epochs. 
The EVPA corrections for all the epochs are shown in 
Table [2] where in column (1) we give the observing code 
of the epoch and list the epochs that were used to anchor 
the D-terms. The epoch of observations is listed in col- 
umn (2), and the reference antenna used in the calibra- 
tion in column (3). The EVPA corrections at 15.4, 12.1, 
8.4 and 8.1GHz are given in columns (4)-(7). Since we 
are using five different anchoring epochs with different 
VLA calibration sources and additionally the UMRAO 
data, the main source of error in our calibration method 
should be the scatter in the measured D-term phases. By 
calculating the standard deviation of the mean for the 
scatter in each right-hand or left-hand phase and taking 
the maximum value over the frequency band as a conser- 
vative error estimate, we determine the absolute EVPA 
calibration errors to be 3°, 2° and 4° at 15, 12 and 8 GHz 
bands, respectively. The total error in the EVPAs is a 
quadrature sum of the calibration error and statistical 
error in the EVPA, with the latter being derived from 
the rms values in Q and U maps. 

The error in the final rotation measure images is highly 
dependent on the error of the EVPA. We have performed 

^ http:/ /www. aoc.nrao.edu/~smyers/calibration/ 
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TABLE 1 

Sources and their rotation measure properties 



lAU Name 


Other name 


z 


Opt. CI. 0app 


Epoch 


Gal 


. RM 


med. RM med. core RM 


med. jet RM 








(c) 




(rad 


m-2) 


(rad ni ^) 


(rad m ^) 


(rad m^^) 


(1) 


(2) 


(3) 


(4) (5) 


(6) 




(7) 


(8) 


(9) 


(10) 


0003-066 


NRAO 005 


0.3467 


B 8.4 


2006-07-07 




4.5 


-20.7 


-34.8 


130.0 


0003+380 


S4 0003+38 


0.229 


Q 


2006-03-09 




-90.0 








0003+380 


S4 0003+38 


0.229 


Q 


2006-12-01 




-90.0 


6053.3 


6053.3 




0007+106 


III Zw 2 


0.0893 


G 1.2 


2006-06-15 




-3.4 


604.2 


604.2 




0010+405 


4C +40.01 


0.256 


Q 


2006-04-05 




-77.8 








0010+405 


4C +40.01 


0.256 


Q 


2006-12-01 




-77.8 








0016+731 


S5 0016+73 


1.781 


Q 8.1 


2006-08-09 




-9.1 


264.9 


264.9 




. — Columns 


are as follows: 


(1) lAU Name (B1950.0); (2) alternate 


name; 


(3) redshift; (4) optical 


classification 


where Q — qua 



arent speed used in viewing angle calculation of Fig. ^ (6) epoch of the RM 



Note. 

BL Lac object, G — active galaxy, and U — unide ntified; (5) Ap 

observation; (7) Galactic RM correction taken from jTaylor et al.| ( |2009[ l; (S) median RM over the source; (9) median RM over the core area; (10) 
median RM over the jet. Table 1 is published in its entirety in the electronic edition of the Journal, a portion is shown here for guidance regarding 
its form and content. 



detailed simulations to verify that our error estimate, 
derived with error propagation from the rms in Q and 
U images, is correct. These simulations are described in 
detail in Appendix [B] where we also give the equations 
used in the error calculation. 

2.3. Image alignment 

During the initial data reduction process the absolute 
coordinate position of the source is lost and the center 
of the image is shifted to the phase center of the map. 
This may not be the same position on the sky for dif- 
ferent frequency bands, and therefore an extra step is 
needed to align the images. This can be done using bright 
components in the optically thin part of the jet, whose 
posi tion should not depend on the observing frequency 



e.g.|Lobanov|1 998MMa rr et al.|200T||Kovalev et al.|2008 

Sokolovsky et al. 2011)^ 'i'his approach works well tor 
knotty jets but is unreliable or impossible to use for faint 
or smooth jets. A solution is to use a 2D cross-correlation 
algorithm to look for the best alignment based on corre- 
lation of t he optically thin parts of the jets at differen t 
bands (e.g. [Walker et al.l[2000{ [Croke k 'Gabuida|[2008| . 

We used both methods whenever possible, and con- 
cluded that the results matched very well when using 
bright opti cally thi n component s. Similarly to e.g., Marr 
et al. ( 2001j and Kovalev et al. ( 2008| ), all the shifts were 
verified by examining the spectral index maps before and 
after the alignment. In shifted maps the spectral index 
gradient along the jet was typically smoother, and any 
optically thin regions apparently upstream of the core 
disappeared. The absolute shifts between 15 GHz and 
other bands varied between Omas and 2.02 mas with a 
median value of 0.11 mas. This is comparable to the pixel 
size of 0.1 mas used in the RM images. The extreme value 
of 2.02 mas is for the source 2134-1-004 between 15 GHz 
and 12 GHz, where a different component is the brightest 
feature in the two maps. This illustrates the importance 
of correct alignment for the data analysis. The small 
median shift, however, shows that in the majority of the 
sources the change is not extreme as is to be expected 
for bright, core-dominated objects. These shifts are de- 
termined in part by the frequency-dep endent core shifts, 
which will be studied in our sample in |Pushkarev| ( j2012j ), 
although other effects can also contribute m some cases, 
such as in 2134-1-004 described above. 

For 35 sources we were not able to find a reliable align- 
ment due to the compactness of the source or the faint- 
ness of a featureless jet. In these cases we aligned the 



images based on the fitted core component position at 
each band. The median shift values for these sources 
were less than 0.03 mas. We used spectral index maps 
to verify that our alignments were reasonable. The spec- 
tral index maps of all the sources will be presented and 
discussed in a separate paper (T. Hovatta et al. 2012 in 
prep.). 

Additionally we did several tests, described in Ap- 
pendix!^ to study the effect of false alignment on spec- 
tral index and rotation measure maps. Based on the 
tests we conclude that even if our image alignment is 
off by 0.15 mas between 15 GHz and any other frequency 
band, it should not affect the results from our rotation 
measure maps, especially as we are not using the edge or 
low signal-to-noise regions to make conclusions about the 
RM structure. We verify that the spectral index map is a 
good indicator of the image alignment because the effect 
of small fake shifts can readily be seen in the structure 
of the spectral index map. 

2.4. Rotation measure maps 

For the calculation of the rotation measure (RM) maps 
we wrote a Perl Data Language (PDL) script that does 
the calculation semi-automatically for our large sample 
of sources. We verified the performance of the script by 
using the RM task manually in AIPS for several sources. 
In our calculations we blanked all the pixels which had 
polarized flux density less than three times the polariza- 
tion error, defined in AppendixiBl at any of the frequency 
bands. Our script chooses the best A^-fit based on a 
criterion and blanks all the pixels where it is not met. We 
calculate the of the fit using the standard formulae 



N 

E 



(O. ~ 



(2) 



where N is the number of data points, Oi are the observed 
data, Ei are the expected data based on the model and 
(Ti is the rneasur ement error of the individual data point 
(e.g. Press|1992 ). Due to the dependence on the errors of 



the data points, blanking of low signal-to-noise regions is 
essential to prevent small values simply due to large 
error bars. In general care must be taken in determining 
the EVPA errors because too small errors will prevent 
good fits while too large errors will result in too small 
values. Our EVPA errors are estimated by adding 
in quadrature an rms error using error propagation from 
Q and U images (see Appendix [b] for details) and an 
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TABLE 2 

EVPA CALIBRATION CORRECTIONS FOR ALL THE EPOCHS IN DEGREES. 



Obscode 


Epoch 


Ref. Ant. 


8.1 GHz 


8.4 GHz 


12.1 GHz 


15.4 GHz'" 


BL137A 


2006 Feb 2 


PT 


-15.9 


-16.7 


-16.1 


-18.8 


BLISTB'^ 


2006 Mar 6 


FT 


-17.7 


-20.0 


-18.0 


-14.5 


BL137C 


2006 Apr 5 


KP 


19.2 


13.1 


22.5 


30.9 


BL137D 


2006 Apr 28 


FD 


12.4 


6.7 


-42.9 


-53.7 


BL137E 


2006 May 24 


FD 


-12.8 


-17.4 


-16.7 


-47.3 


BL137F'' 


2006 Jun 15 


FD 


-42.2 


-47.3 


-47.1 


-49.2 


BL137C 


2006 Jul 7 


FD 


-46.3 


-47.9 


-47.9 


-49.1 


BL137H 


2006 Aug 9 


FD 


-45.0 


-47.8 


-47.0 


-48.4 


BL137P 


2006 Sop 6 


PT 


-10.0 


-10.0 


-14.3 


-14.9 


BL137J 


2006 Oct 6 


FD 


-45.3 


-46.8 


-45.6 


-48.3 


BL137K'' 


2006 Nov 10 


FD 


-44.0 


-45.0 


-45.0 


-46.6 


BL137L 


2006 Dec 1 


FD 


-44.2 


-47.1 


-47.4 


-50.5 



^Epoch used to anchor the D-term caUbration 
'^Cahbrated as part of the MOJAVE project 

absolute calibration error defined in Sect. 12.21 As we 
are fitting a two-parameter model to four data points we 
have two degrees of freedom and from a distribution 
the corresponding 95% confidence limit is < 5.99. 

The EVPA is ambiguous for changes of 180° and in 
the calculation of the RM we need to solve for these 
nTT-wraps. We first assumed that there are no nvr-wraps 
between our frequency bands and calculated the RM fit. 
If the of the fit met our criterion, we accepted the 
RM value without any wraps. If the x^ criterion is not 
met, we solved for all possible nvr-wraps up to 3.3 x 10"' 
rad m~^ and chose the fit with the smallest wrap meet- 
ing our x^ criterion. The upper limit was primarily intro- 
duced to keep the computing time reasonable but also be- 



cause based on earlier studies (e.g., Zavala & Taylor 2003 
2004 1 , we did not expect to resolve RMs larger than this 



with our frequency setup. If none of the wraps resulted 
in acceptable fits, we blanked the pixel. By blanking the 
poor A^-fit regions, we prevent interpretations based on 
noisy data and identify regions with non-A^-law behav- 
ior. 

The error of the RM is calculated from the variance- 
covariance matrix of the least squares fit in each pixel. 
Our typical errors range between 70 and 150 rad m~^ 
depending on the signal-to-noise ratio of the total in- 
tensity in the jet, thus in the fainter jet edges the RM 
errors are larger. We verified that our error estimates are 
correct by performing detailed simulations described in 
Appendix [C] 

In order to study the distribution of the intrinsic, 
redshift-corrected, RM values, the Galactic Faraday ro- 
tation contribution must be taken into account. We used 



the averaged Galactic RM image of Taylor et al. ( 2009 1 
and subtracted the value at the source location trorn 
each map. We list the values used for each source in 
Table [TJ In the majority of sources the Galactic Fara- 
day rotation is very small (the median absolute value for 
the sample is 12.3 rad m~^) but there are 14 sources 
for which the absolute value is more than 70 rad m~^, 
thereby exceeding our minimum error in the RM values. 
The largest absolute Galactic Faraday rotation values are 
observed for 2021-1-317 (-173 rad m^^) and 2200+420 
(—156 rad m^^). For the majority of sou rces in our 
sample the values from Taylor et al. (|2009[) agree very 
well with previously pubTTshed values (iRudnick fc Jones 
19831 |Rusk||1988| |Wrobel||1993l |Fushkarev||200l|). How- 



ever, we note that smce we are usmg an averaged image. 



some Galactic RM values may be underestimated be- 
cause small regions of high Galactic RM get smoothed 
out (e.g., 0235-M64, 1749-1-096, 1803-^784, 2200-^420). 

3. RESULTS 

RM maps are shown in Fig. [TJfor the 159 cases where 
we detect enough polarization to get a RM value for at 
least a few pixels. We show the RM values in color scale 
overlaid on the 15 GHz total intensity contours and ex- 
amples of the fits in two locations of the jet, chosen 
to be at the polarization peaks of the map. These lo- 
cations are typically at least one beam width apart. In 
sources where clear polarization peaks were not seen, we 
chose the location to be in the middle of the RM region. 
Additionally we show the error of the RM in color and 
the intrinsic, RM corrected, 15 GHz polarization vectors 
overlaid on the 15 GHz polarization contours. All the RM 
maps in Fig. [T] and later in the paper are corrected for 
Galactic Faraday rotation. In some cases, there appear 
to be pixels with very high RM values of over ± 2 x lO'' 
rad m^^. In most of the sources these coincide with edge 
pixels and/or regions of complex polarization structure. 
Our simulations of the RM error (see Appendix[C|, show 
that it is possible to have spurious high-RM pixels in the 
maps purely due to random noise in the polarization im- 
ages. These were always more than ± 2 x 10^ rad m~^ 
and in our real maps the high-RM regions resembled the 
simulated maps very well. Therefore we have blanked 
these extreme values in Fig. [l] 

3.1. Extreme RM values 

In some sources we see very high RM regions around 
the core where we may expect more F araday rotating ma- 
terial and stronger magnetic fields. |Udomprasert et al.| 
( |1997| report an intrinsic RM (RMi„t ) of 4 x lO'' rad m"^ 
in the high-redshift quasar OQ 172 (z = 3.52). In the 
observed frame, defined as RMobs = (1 + •z)^^RMint, 
this corresponds to a RMobs of ~ 2000 rad m~^ in the 
core. If the intrinsic value is correct, we might expect 
to observe extremely high RMs in some nearby objects. 
Attridge et~al] ( |2005| ) report a difference of ~ 3.2 x 10"* 
rad m~^ between two components in the core of 3C 273 
in observations at 43 and 86 GHz. It is, however, very 
difficult at our observing frequencies with much less res- 
olution to distinguish true extremely high RMobs val- 
ues from the spurious ones due to noise and blending 
of components. For example, in 3C 273 we observe ex- 
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Fig. 1. — Rotation measure maps of all the sources. RM map (top left) with A^-fits at one or two locations in the jet (top right). Error 
map (bottom left) and map of rotated EVPA values (bottom right). Figures 1.1-1.159 are available in color in the electronic edition of the 
Journal. 



treme RMobs values of ~ 2.3 x lO"*^ rad m"^ around the 
core in the March epoch. In our June epoch, we do 
not detect these high values but instead see values of 
~ —2.9 X 10* rad around the same region. Similar 
behavior is seen in the cores of 3C 279 and 3C 454.3. In 
2200+420 we do not find good A^-fits in the core in our 
April epoch but in November we detect extreme RMobs 
values of ~ +2.9 x 10* rad m~^, never seen before in this 



source (Mutel et alj 2005 



ing observations of O'Sul 



1[200 



Zavala fc Taylor 2003| includ- 
ivan & Gabuzda' ( 20U9) in July 



2006. Other sources where we detect extreme RMobs 
values in larger areas near the core include 0149+218, 
0420-014, 0605-085, 1038+064 and 2145+045. Out of 
these, 0420-014 was observed by [Zavala fc Taylor] ( [2003 1 
who do not see these extreme values although they too 
do not find a good fit at the core compone i it pos ition. 
0605-085 was observed by [Zavala fc Taylor[ pOOi] ) who 
do not detect any extreme RM values and lind the RM 
in the core to follow the A^-law. In 2145+067 the polar- 



ization structure at 15 GHz is extremely complex with 
four separate components seen within the innermost jet, 
while at 8 GHz only one component is seen. Therefore 
we do not believe that these extreme observed core RMs 
are real in our maps, but instead are due to multiple po- 
larized components blending in the finite beam or due 
to different opacity properties at the frequency bands. 
This is further supported by continuing MOJAVE ob- 
servations of 2200+420 which show a ne w component 
emerging from the core in February 20071^ Algaba et al 



(20111 observe high RMs in the cores ofseveral sources 



in their study of eight sources between 12 and 43 GHz. 
They observe a RMobs of 2.2 x 10* rad m~^ in the core of 
1633+382, for which we observe a small RMobs of —244 
rad m^^. Their result is based on a large difference in the 
EVPAs of the 22 and 24 GHz observations which require 
a large RM. It is possible that they are able to resolve 

* http: / /www. physics. purdue.edu/astro/MOJAVE/sourcepages/2200+420.shi 
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structures not seen in our maps due to their higher reso- 
lution. In some sources they also find that they need to 
divide the frequency range into high- and low-frequency 
parts to obtain acceptable A^-fits, which is further indi- 
cation of different frequencies probing different regions 
in the core and multiple components blending within the 
beam in the lower frequency maps. 

The blending of components in the core region can also 
affect our A^-fits so that the criterion is not met and 
no RM values are shown in the maps. Another cause for 
this could be internal Faraday rotation, which could play 
a significant role in the AGN core regions. We also see 
non-A^ patches in the jets of some sources, sometimes due 
to the faintness of the jet emission but at other times also 
due to depolarization of the lower frequencies, a sign of 
internal Faraday rotation. The effects of internal Fara- 
day rotation and other depolarization mechanisms are 
discussed in more detail in Section l4?2l 



3.2. Median RM distribution 

We were able to determine the median RMobs for 159 
maps, which are shown in the top panel of Fig. [2] and 
given in TablejlJ where column (1) gives the B1950-name 
of the source and column (2) an alternative alias name. 
The redshift and the optical classification of the source 
are listed in columns (3) and (4). Apparent speed used 
for calculation of the de-projected distance in Sect. 4.1 is 
given in column (5) and the observing epoch is listed in 
column (6). The value us ed for Galactic Faraday rota- 
tion correction, taken from Taylor et al.,(2009j is listed in 
column (7). The median KMobs value, taken as the me- 
dian of all the pixels in the source where RM is detected 
and not blanked, is listed in column (8). Columns (9) 
and (10) give the median RM over the core and the jet 
regions, respectively. We calculate the median instead 
of the average to lessen the effect of individual, possibly 
spurious, high-RM values. The vast majority of sources 
have a median RMobs of less than 1000 rad m~^, but the 
distribution has a tail to RMobs values of 6500 rad m~^. 
The highest value shown in the plot, 6457 rad m^^, is 
for the source 2008—159, which only shows RM values 
in a small region of less than half the beam size. At the 
redshift 1.18 of the source, this would result in an ex- 
tremely high intrinsic RMint of over 3 x lO** rad m~^ in 
the source frame. As the region over which we detect the 
high RMobs value is so small and does not coincide with 
any total intensity component locations, it is difficult to 
say if this is a true RM of the source or due to blending 
of multiple components within the core region. 

3.3. Core vs jet distributions 

To study the difference between core and jet RM val- 
ues we 1) divided the source into core and jet regions by 
defining the core region to be everything within a beam 
width from the center of the 15 GHz core component po- 
sition and jet region to be everythin g else and 2) took the 
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15 GHz modelfit components (see 
for details on the modclfitting) and divided the source 
into core and jet components. In the first approach the 
division is determined by a line perpendicular to the jet 
direction at one beam width away from the core. Median 
values over the pixels within the regions were calculated 
and are given in Table ^columns (9) and (10) and shown 
in the bottom two paiiels of Figj2[ 



^^^1 Un— identified - 
^^<^<^ Galaxy 

BL Lac object 




1000 2000 3000 4000 

Median iRM^j^gl over the whole source (rad/m^) 



Un— identified - 
Galaxy 

BL Lac object 
Quasar 




1000 2000 3000 4000 

Median |RM^ | in core area (rad/m^) 



Un— identified 
Galaxy 

BL Lac object 
Quasar 



Median iRM^^^gl in jet area (rad/m^) 

Fig. 2. — Distribution of median |RMobsl over the whole source 
(top), over the core region (middle), and over the jet region (bot- 
tom). Quasars are shown in white, BL Lacs in hatched, galaxies 
in cross-hatched and un-identified sources in black. These plot ex- 
cludes quasars 0003-1-380 and 2008-159 which have RM values over 
6000 rad m~^. 



In the second approach, we calculated the average RM 
over the 9 contingent pixels around the component po- 
sition to avoid basing conclusions on single pixel values. 
This corresponds to 10-30% of the restoring beam width 
depending on the declination of the source. The compo- 
nent locations and their RMobs values are given in Table 
[Sj The I.D. number of the component is listed in column 
(2) where indicates a core component. Columns (4) 
and (5) give the component distance and position angle 
from the phase center of the map. The RMobs and its 
error are given in column (6). Because the pixels are not 
independent (i.e., they cover a region smaller than the 
FWHM of the restoring beam), it is not straightforward 
to estimate the error on the average carried out over 9 
pixels. We define the error as the average of RM errors 
in the 9 individual pixels; this approach is conservative, 
and may overestimate the true error somewhat. Column 
(7) indicates whether a jet component is isolated (see be- 
low). Most of the sources in the MOJAVE sample are 
core dominated, with a bright compact core that is opti- 
cally thick at centimeter wavelengths and a fainter jet. In 
most of our sources, we identify the core as a bright, sta- 
tionary feature in the jet, typically at one extreme end of 
the jet. In a few sources (especially with two-sided jets) 
the identifica tion is not as simple and these are discussed 
separately in Lister et al. (2009b). 

The distribution ot the RMobs values in the compo- 
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TABLE 3 

MODELFIT COMPONENTS AT 15.3 GHz AND THEIR, RMobs VALUES 



Source 


I.D. 


Epoch 


r 


P.A. 


RM 


Isolated 








(mas) 


(deg) 


(rad m~'^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0003-066 





2006-07-07 


0.71 


-168.6 


-78 ± 73 




0003-066 


1 


2006-07-07 


0.66 


-71.6 


-88 ± 75 




0003-066 


4 


2006-07-07 


6.86 


-81.1 


113 ± 98 


Y 


0003-066 


5 


2006-07-07 


0.09 


2.9 


-42 ± 72 




0003-066 


6 


2006-07-07 


1.28 


-102.7 


-14 ± 77 




0007+106 


1 


2006-06-15 


0.39 


-66.6 


627 ± 130 




0016+731 





2006-08-09 


0.01 


-52.2 


278 ± 74 




0016+731 


2 


2006-08-09 


0.20 


129.2 


266 ± 74 





Note. — Columns are as follows: (1) lAU Name (B1950.0); (2) I.D. of the component (0 = core); (3) Observing epoch; (4) Component 
distance from the phase center of the I map; (5) Position angle of the component from the phase center; (6) Component RM; (7) Flag for 
isolated jet components. Table 3 is published in its entirety in the electronic edition of the Journal; a portion is shown here for guidance 
regarding its form and content. 



nents is shown in Fig. 3] We were able to determine 
the core component RM in 104 maps (101 sources) and 
the jet RM in 324 components (121 sources). From the 
distributions it is clear that the core component values 
have a tail to higher RMobs values, but there are also 
some jet components with high RMobs values. In most 
of those cases the jet component is within 1 mas of the 
core component and often still in the optically thick or 
self-absorbed region of the jet. To distinguish the jet 
components which are away from the bright core region, 
for each component we calculated the combined contri- 
bution of all the other jet components in the map at the 
component's peak intensity position. If this sum was less 
than 30% of the component's total intensity, we consid- 
ered the component to be isolated. In this way we deter- 
mined that the polarization and RM of the component 
were not affected by nearby bright components. Out of 
all the jet components, 36 in 24 sources are listed as iso- 
lated. The distribution of RMobs values of these isolated 
components is plotted in the bottom panel of Fig. [Sj As 
can be seen, none of the RMobs values greater than 700 
rad m~^ are isolated. The median component RMobs 
in the whole sample is 171 rad m^^ for the cores, 125 
rad m~^ in all the jet components and 104 rad m~^ in 
the isolated jet components. Accor ding to an A nderson- 



Press|1992 '), which is 
ran the Kolmogorov- 



Darling (A-D) two sample test (e.g., - 
more sensitive to distribution tails tl 
Smirnov test, the core and jet components have less than 
a 1% probability of coming from the same parent popu- 
lation. Comparison of core and isolated jet components 
gives a probability of less than 2% due to the smaller 
number of jet components. In all our tests we consider 
the result significant if the probability is less than 5%. 

A similar trend is seen when comparing the core and jet 
regions m Fig. The median RM for the core regions 
is 187 rad m~ and for the jets is 102 rad m~^, a re- 
sult very similar to that for the isolated jet components. 
According to the A-D test the probability for these dis- 
tributions to come from the same parent population is 
less than 0.001%. 

3.4. Optical Subclasses 

Our sample can also be divided into subclasses based 
on the optical classification of the source. These are also 
shown in Figs. [2] and [3j The number of galaxies and 
optically un-identified sources in our sample is so small 
that they cannot be included in any statistical compar- 



1 1 1 1 1 1 1 1 
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Fig. 3. — Distributions of IRMo^sl ™ the model fit core (top) 
and jet (middle) components. Jet components which have less 
than 30% total intensity contribution from other components in 
the source are plotted in the bottom panel (see text for details). 
Quasars are shown in white, BL Lacs in hatched, galaxies in cross- 
hatched and un-identified sources in black. 



isons. The quasars and BL Lac objects, however, can 
be compared, and it is clear that the high-RM tail in the 
distributions consists mainly of quasars. The median ab- 
solute RMobs value in quasars is 144 rad m^^ and in the 
BL Lacs 79 rad m~^. 

If we look at the core and jet components individually, 
there is less than 0.1% probability that the jet compo- 
nents of quasars and BL Lacs are drawn from the same 
population. The median jet RMobs for quasars is 141 
rad YcT^ while for the BL Lacs it is 71 rad m~^. How- 
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ever, this difference is affected by tlie components witliin 
one beam width of the core in quasars because the me- 
dian values in the jet regions of quasars is 116 rad m~^ 
and in BL Lacs it is 76 rad m~^. According to the A-D 
test we cannot reject the null hypothesis that the distri- 
butions come from the same parent population. 

In the cores differing results are also obtained when 
core components and regions are compared. In quasars 
the median core component value is 183 rad m~^ com- 
pared to 134 rad m"^ in BL Lacs, and we cannot reject 
the null hypothesis that they come from the same parent 
population. In the core regions, however, the median for 
quasars is 200 rad m~^ and for BL Lacs it is 105 rad m^^, 
and there is less than 1% probability that these come 
from the same distribution. 

The se results can be compared to [Zavala fc Taylor] 
(2004), who saw a difference in the core RMs for quasars 
and BL Lacs but not in the jet values. Our results on 
the jet RMs agree if we look at the jet regions which 
are not contaminated by components near the core. We 
cannot verify if there is a difference between the cores 
in quasars and BL Lacs because of the differing results 
depending on if we look at the core components or core 
regions. However, the higher median values observed in 
quasars than in BL Lacs are in accordance with the stan- 
dard models in which BL Lac objects have less material 
around them, resulting in dimmer and narrower emission 
lines than in quasars. 

When examining the intrinsic, redshift corrected, 
RMjiit values in the components the difference between 
BL Lacs and quasars is more significant. This is mainly 
due to the BL Lacs in our sample having smaller red- 
shifts than the quasars (median redshift 0.31 vs. 1.12), 
which is enhanced in the correction defined as RMjnt = 
RMobs(l + z)'^, where z is the redshift. 

The absolute RMj^t values in the cores of quasars range 
from 4.8 to 6436 rad m~^ with a median of 798 rad m~^. 
In BL Lacs the range of the core RMint values is from 13 
to 3873 rad m~^, similar to the quasars, but the median 
is significantly smaller (274 rad m~^) and according to an 
A-D test the probability for the two distributions to come 
from the same parent population is less than 3%. The 
median RMint value for the jet components of BL Lacs 
(range from 0.8 to 1937 rad m~^ with a median of 148 
rad m^^ ) is significantly smaller than that of the cores 
even though the range is similar. In quasars they range 
from 1 to 8975 rad m^^ with a median of 563 rad m^^ but 
the probability to reject the null hypothesis in the case of 
the intrinsic jet and core components of quasars is only 
6.8% and therefore not significant. Similarly to the case 
of core components, the difference between quasars and 
BL Lacs is significant also in the intrinsic jet components. 

4. DISCUSSION 

One of the main scientific motivations for the multifre- 
quency survey of the MOJAVE sources was to determine 
the effects of Faraday rotation on the observed polar- 
ization structure of the sources at 15 GHz. Based on 



the first epoch MOJAVE data, [Lister fc Homan ( [2005 1 
showed that in BL Lac objects the distribution of Evi^As 
with respect to the local jet direction appears bimodal. 
The effect of Faraday rotation was not taken into ac- 
count and therefore these results could be affected by 
sources with high RM values. The RM distribution of 



Fig. [2] shows that in over 80% of our sources the RMobs 
values are less than 400 rad m~^, which will rotate the 
15 GHz electric vectors by about 10 de grees. This means 
that t he results of a large sample of [Lister fc Homan] 
(2005) should approximately reflect the true distribu- 
tion at 15 GHz. However, when studying some individual 
sources, the Faraday rotation must be taken into account 
as a rotation measure of 2000 rad m~^ (seen in the me- 
dian RM distributions for a few individual sources) can 
rotate the 15 GHz EVPAs by 40 degrees. For example, 
in 0429-f415 we detect RMobs of 1900 rad m'^ in sev- 
eral jet components 40 mas from the core, similarly to 
Mantovani et al.[ ( [2010[ ). In 1101-1-384 and 1725-1-044 we 
detect core component RMobs as high as 3800 rad m^^, 
although it must be noted that in these two sources we 
detect RM only in a very small region around the core 
and therefore cannot be sure if it is true RM or due to 
blending of multiple polarized components within the fi- 
nite beam. 

4.1. Distance dependence 

In Fig. [4] (top panel) we show the RMobs versus pro- 
jected distance from the core for all the modelfit jet 
components. The dependence for the total sample is 
not very clear although according to a non-p arametric 



Kendall's tau correlation test (e.g., Press|T992 ) there is a 
significant negative correlation (t = —0.13, p=0. 00058). 
When quasars and BL Lacs are studied separately it 
can be seen that the correlation in quasars is stronger 
(r = -0.17, p=6.3 X 10"^) while for the BL Lacs 
alone r = —0.21, p=0.02. The picture is, however, 
more complicated as the true distance from the core 
depends on the viewing angle of the source. We have 
therefore de-projected the distances using viewing angles 
e = tan-i[(2/3app)/(/3^pp + Dyar^ - 1)] determi ned with 
Doppler factors Dvar from Hovatta et al.' ('2009|) and ap- 
parent speeds 



from Lister et al. (2009b). Some of 



the speeds have been updated since iListcr ct aL]j2009b), 
and all the speeds used are tabulated in Table [1] Both 
values were available for 138 components in quasars, 47 
components in BL Lacs and 4 jet components in galaxies. 
The RMobs against the de-projected distance is shown 
in Fig. [4] (bottom panel). The negative correlation in 
quasars remains significant despite the smaller number of 
sources (r — —0.23, p=4.8 x 10"^) while in BL Lacs the 
correlation vanishes (r ~ —0.01, p=:0.92). However, we 
do not detect jet components as far away from the core in 
BL Lacs as in quasars which could affect the correlation. 
The correlation seen in quasars supports the results from 
the simple core and jet component comparison in Sect. 
[3.3| showing that the amount of Faraday rotating mate- 
rial diminishes as a function of distance from the core. 

The above approach is a simplification of a more com- 
plex behavior in the RM values as a function of distance 
and even if the general trend shows a decline in RM along 
the jet, individual sources may deviate from this trend 
and show complex structures. Another way to study the 
distance dependence is to calculate the RM values along 
the total intensity ridge line of the jet. Two examples 
are shown in Fig. [5] where the top panel shows the RM 
along the ridge line in 0430-1-052 where the core is depo- 
larized (Fig. 1.28) but further along the jet RM declines 
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Fig. 4. — IRMobsl against the projected distance of the component from the core (top panel) and against the de-projected distance 
taking the viewing angle of the source into account (bottom panel). Left panels show quasars (crosses), and right panels BL Lac objects 
(asterisks), and galaxies (triangles). 
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Fig. 5. — RMobs along the total intensity ridge line in 0430+052 
(top) and 2230+114 (bottom). The average FWHM beam size 
along the jet is shown in both plots as a scale bar. 



very sharply, in accordance with the simple scenario. In 
2230 +114 (Fig. 1.151) a more complex structure can 
be seen along the jet, with the RM changing sign along 
the jet. This example shows that with better resolution 
along the jet, the situation may not be as simple as a lin- 
ear dependence along the jet but more complex regions 
are seen. The recent sensitivity upgrades of the VLBA 
(e.g., higher bit rate observations) will allow us to detect 
more polarization further down from the core and help 
to study this in more sources. 

4.2. Faraday depolarization 

Faraday rotation can cause different amounts of depo- 
larization depending on the nature of the rotating screen 
and if it is internal or external to the jet ( Burnj [l966). 
By applying a very simple Burn-type internal I'araday 
depolarization model to the core components of 40 AGN, 
Zavala & Taylor ( 2004 ) concluded that internal Faraday 
rotation alone cannot explain the steep decline in frac- 
tional polarization as the magnitude of the rotation mea- 
sure increases. The equations they used are valid only 
in the optically thin regime and therefore not applicable 
to the core regions of AGN at 15 GHz. We explore the 
viability of possible models by directly fitting individual 
isolated jet components for depolarization and compar- 
ing the results to our observed RM values. 

For internal depolarization assuming a uniform mag- 
netic field and the optically thin regime we have 



mobs(%) 



sin(2A2RM) 



2A2RM 



(3) 
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where niobs is the observed fractional polarization, mmax 
is the maximum fractional polarization in the specific 
magnetic field configuration rl RM is the observed rota- 



tion measure and A is t he observing wavelength (Burn 
1966 Homan et al.|2009 ). In case of external depolariza- 
tion 



mobs(%) = 



(4) 



where a is the standard deviation of the RM fiuctu ations 
and the rest of the parameters are as in Eq. [S] ( [Burn 



19661. Here we assume that the component is optically 



thin and homogeneous (not a combination of multiple 
components), and also that the angular scale of RM vari- 
ations is much less than the angular resolution of our 
observations. The functional forms of depolarization in 
Eqs. [3]and|4] are similar over the range of |RMobs| we ob- 
serve in the isolated jet components (up to 800 rad m^^) 

and both follow the functional form m = moe where b 
is — 2RM^ in the case of internal depolarization and — 2cr^ 
in the case of external depolarization. We can linearize 
the formula to fit In m = In mo -I- bX'^ to our observations. 
This way we will get an estimate of total depolarization 
b from the slope of the fit and the intercept, lump, gives 
the maximum polarization for that component. We use 
the isolated jet components only to ensure that we are 
looking at homogeneous components in the optically thin 
part of the jet. The polarization values for the isolated 
components are given in Table H] where the RM is given 
in column (6), columns (7) - (10) show the fractional po- 
larization and its error at the different frequency bands, 
and column (11) the value b. In Fig.p^we show examples 
of the fits and in the top panel of Fig. [7] we show the 
square root of |6| with its sign preservedto distinguish 
depolarization from inverse depolarization where polar- 
ization is higher at 8 GHz than at 15 GHz. In 16 out of 61 
components our simple model does not fit the data well 
which may be an indication of a more complex behavior 
than the simple exponential model can explain. These 
components are clearly marked in Fig. [7] and Table |4j 

The solid line in the top panel of Fig. [7] gives the ex- 
pected amount of internal Faraday depolarization for a 
given RM. If internal Faraday depolarization alone would 
account for our observations, most of the data points 
should fall on this line which is not the case. The ma- 
jority of the depolarized components fall below the line 
indicating that the fractional polarization in these com- 
ponents falls faster than expected for internal Faraday 
rotation. This is clearer in the middle panel of Fig. |7] 
which shows only the components that are more than 2a 
away from 6 = 0. The same appears to be true for exter- 
nal depolarization when we assume that the dispersion 
in the RM values is proportional to the observed RM, 
i.e. (T^ — RM^ (dashed line). Most of the components 
fall below this line, indicating that the dispersion in the 
Faraday screen is larger than the mean RM produced 
by that screen. This suggests that the Faraday screen 
is dominated by random RM fluctuations between inde- 
pendent lines of sight. For a random Faraday screen the 
observed average RM will approximately follow a rela- 
tion RM oc cr/\/N where N is the number of lines of sight. 

^ Note that mm ay ~ 70% for a pure uniform magnetic field 
| |Pacholczyk||1970] l; however, we do not assume a value for mmax 
in our analysis and include it as a free parameter, mo, in our fits. 



The dotted line in the top panel of Fig. [7] uses N = 10 
for the calculation of a in Eq. [4j where we assume the 
angular scale of the RM dispersion to be much smaller 
than the beam so that = lORM^. This line fits our 
data much better than assuming ~ KM^ , but there is 
still a large scatter about the line. The line produced by 
Eq. |4] assumes the scale of Faraday rotating cells to be 
much smaller than the beam size. This may not be true 
for high-an gular-resolutio n observations such as these by 
the VLB A (Tribble 1991 1. In order to take the number of 
lines of sight correctly into account, we directly simulate 
the expected depolarization and RM for a variety of cr, 
RM and N combinations. 

In the simulations, we define values in the range of 
to 1.6 X 10^'^m^ and initialize each frequency to have 
70% fractional polarization and 0° EVPA.FJWe pass this 
initial polarization through N individual lines of sight. 
Each individual line of sight is simulated by adding an 
average RM, which is the same for all lines of sight, and 
a random contribution drawn from the Gaussian distri- 
bution of variance cr^. For each wavelength, we sum the 
contribution of N lines of sight drawn in a similar man- 
ner to obtain an average p and EVPA value. We treat 
these average values as our real observations and fit the 
p values for b and EVPA values for RM. We repeated 
the simulation 1000 times to obtain a range of b and RM 
values to compare with our observations. In Fig. [7] bot- 
tom panel we plot several of our simulations and in the 
top panel we show our best fit case overlaid with the real 
data points. 

The blue dots in the bottom panel of Fig. [7] are from 
several simulations with a large number of lines of sight 
(N=1000) where we have set average RM — a (using val- 
ues 100, 200, 300, 400, 500, 600, and 700 for RM and a), 
and these simulations are plotted with the dashed line 
produced by Eq. |4j As can be seen, these large N simu- 
lations follow the expected curve very well with increased 
scatter for the large sigma values. As described above, 
our data in the top panel of Fig. [7] largely fall below this 
line, indicating that the Faraday screens may be random 
screens with a small number of lines of sight. To test this 
possibility, we set the average RM applied to all lines of 
sight to be in our simulation, the number of lines of 
sight N = 10, and either a = 200 rad m~^ (black dots) or 
a = 300 rad m~^ (red dots). The red dots cover almost 
the same region as our data while the black dots pro- 
duce too little depolarization. As one might expect for 
a random- walk style Faraday screen, ct/\/N for the red 
dots is 95 rad m~^ which is in close agreement with the 
observed median |RM| — 104 rad m~^ for our sample of 
isolated components. Therefore we conclude that most of 
our observations can be explained with a completely ran- 
dom external foreground screen viewed through a small 
number of lines of sight. This implies that the linear size 
of the RM cells may not be too much smaller than our 
beam size. We plot the red dots of the bottom panel of 
Fig. [7] in the top panel to show the good correspondence 
with our data. 

Additional complications to note are that if the depo- 
larization is much higher, we do not detect enough po- 

^ Note that the initial values chosen here do not affect our results 
as we are interested only in the change in p and EVPA due to 
combination of multiple lines of sight within our beam. 
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Fig. 6. — Examples of fits of depolarization curves to isolated components, a) No apparent depolarization, b) Ambiguous behavior, c) 
and d) Depolarization, e) and f) inverse depolarization. See text for details. Upper limits are shown in arrows. Note the scale of the 
vertical axis which is not linear but in natural logarithm. Fits for all isolated components are available in the online edition of the Journal. 



larization at 8 GHz to calculate the RM or the internal 
rotation causes non-A^-law behavior and we do not get 
good enough fit in our RM calculation. To study this 
further, we have included in Table |4] isolated jet com- 
ponents for which we detect fractional polarization at 
some of our four frequency bands but not necessarily all 
(15 components), and ten components for which we de- 
tect fractional polarization but no RM (a sign of non-A^ 
behavior). In the calculation of the slope b, we have as- 
sumed the upper limits to be detections. This way we 
get a lower limit estimate for the depolarization. 

Based on the polarization behavior in Fig. [7] we can 
divide all our fits into four categories. Constant polar- 
ization over the frequency range is seen in 13 out of 60 
components, as the example case in Fig. [6^. These com- 
ponents are within error bars of zero b in Fig. [71 In 18 
components, the fractional polarization did not follow a 
linear trend but was changing randomly between the fre- 
quencies, as seen in Fig. ^p. In several of these the slope 
b is consistent with zero within the error bars. Depolar- 
ization is seen in 20 of the components, and two examples 
are given in Figs. [6]: and [6ji. In 1828-1-487 we were not 
able to calculate a RM value because we only detect an 
upper limit in linear polarization at the 8.1 GHz band. 



From the slope of the fit we can estimate that the amount 
of internal Faraday rotation required to cause such de- 
polarization would be 970 rad m~^, higher than what 
we observe in any of the isolated components. There are 
four additional components which show slopes steeper 
than the typical range in Fig. [TI and the depolarization 
in these sources could be produced by internal Faraday 
depolarization. 

Additionally, we see nine components with inverse de- 
polarization structure, where the fractional polarization 
at 8 GHz is higher than at 15 GHz. In only five of these 
we detect RM as well and these are the most significant 
points above zero 6 in Fig. [7| The other inverse depolar- 
ization components above zero are in the category where 
a linear fit did not describe the fractional polarization 
behavior well and the slope 6 is not a good indicator of 
the depolarization. The nine components each show a 
significant rise in the fractional polarization as shown in 
Figs. [6fe and [6f. This is unexpected and cannot be eas- 
ily explained with any standard external depolarization 
models. Interestingly, seven of the nine components (and 
all the five for which we have RM value) are in 3C 273 
and in 3C 454.3, both of whic h sh ow transverse RM gra- 
dients in their jets (see Sect. 4.3). The other two are in 
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Fig. 7. — Top: Fitted slope b from Fig. ^against the |RMobsl in 
the isolated jet components. Crosses (black in the online version) 
have good fits with < 6, open circles (blue) have x'^ > 6 and 
may not be well described by our simple model. Grey dots (red 
in the online version) are our best-fit simulation case from the 
bottom panel. The solid line shows the expected amount of internal 
Faraday depolarization and dashed and dotted lines the expected 
amount of external Faraday depolarization for different parameters 
(see text for details). Middle: Same as top panel but showing 
only components that are more than 2u from fe = 0. Bottom: 
Simulations for external depolarization over varying average RM, 
number of lines of sight N and dispersion c in the RM. The dashed 
line shows the expected amount of external Faraday depolarization 
when average RM = c. Dark gray (blue in the online version) dots 
show the same from several simulations using N=1000. Black dots 
have average RM = 0, N = 10, a = 200, light gray (red) dots have 
average RM = 0, N = 10, cr = 300. 



1458+718 where the slope is still within 2a from zero, 
and in 1514—241 where the fractional polarization rises 
from an 8.6% upper limit at 15 GHz to 24% at 8 GHz. 
In this source we do not detect any RM values. Internal 
Faraday rotation together with helical or loosely tangled 
random magnetic field configurations could possibly ex- 
plain the observed inverse depolarization and this model 
is investigated in detail by Homan (2012). 



4.3. Transverse RM gradients 

If AGN jets are launched from a rotating black hole 
or accretion disk, it could be expected that the mag- 
netic field around the jet has an ordered toroidal compo- 



nent (e.g., McKinney & Narayan 2007). A signature of 
such a toroidal component (often mterpreted as a com- 
ponent of a helical field) would be a rotation measure 
gradient transverse to the jet flow direction as t he line 
of-sight magnetic field changes its direction (e.g. 



Bland- 



[ford,1993) . In this case, the gradient should be seen m 
multiple locations of the jet, which distinguishes it from 
isolated local gradients that arise from changes in the 
density of the Faraday rotating material. The detection 
of such gradients is challenging due to the limited num- 
ber of bright sources w ith polarized, well-resolved jets 
( Taylor &: Zavala]|2010 ). Furthermore, the jet structures 
can be very complex, and it is likely that both kind of 
gradients exist in th e same sources, as i n the case of the 



radio galaxy 3C 120 ( Gomez et al.|20lT |. Therefore even 
if a transverse RM gradient is observed, it does not auto- 
matically indicate the presence of a helical magnetic field, 
and detailed m odeli ng is needed to probe its nature. 

In Appendix |C.2[ we perform simulations to investi- 
gate how large spurious transverse gradients can arise 
due to image noise and finite restoring beam size. Based 
on our simulations we conclude that the convolved jet 
should be at least 1.5 beams wide (but preferably more 
than 2) in polarization along the direction of the gradi- 
ent and that a gradient should exceed the 3(t level to be 
considered significant. We define a as the largest RM 
error at the edge of the jet when the systematic error 
due to absolute EVPA calibration, ^ 60 rad m~^, is first 
removed in quadrature. The significance of a gradient is 
then simply the total change in RM divided by the a. 
These criteria are s i milar to th e ones described by|Brod- l 
erick fc McKinneyj ( |2010[ ) and [Taylor fc Zavalaf ( |20To| ; 



although our simulations indicate a minimu m transverse 
width of 1 . 5-2 rat her than 3 beamwidths. IBroderick fcl 
McKinney (2010) show that due to the complexity of 
AGIN cores, gradients within one beam width of the core 
may be unreliable at our resolution, and ther efore we 
have not cons i dered these regions in our study. |Murphy| 



& Gabuzda (2011) argue that a RM gradient due to 
helical magnetic field is significant even when the jet is 
not resolved, based on simulations where they convolve 
a simulated gradient with different beam sizes. However, 
their simulation does not take into account that a spu- 
rious gradient can arise due to noise in the data, which 
we show to be a major effect on VLBA observations of 
unresolved jets. 

Following the above guidelines, we examined all our 
RM maps in detail. Our observations show a clear gra- 
dient across the jet of 3C 273 (F i g. confi rming the 



observations o f Asada et al. ( 2002 



20D8a) and Zavala &; 



Taylor (|2005|). I'he gradient is detected above the '6a 
level, and the jet is nearly 3 beams wide along the gradi- 
ent direction. For the first time the RM is seen to change 
sign over the gradient, which is a further indication of a 
helical field. We believe we are seeing this now due to a 
different part of the jet being illuminated in the earlier 



observ ations, similarly as seen in 3C 120 by Gomez et al. 
(2011 ). In fact, if we compare the mean jet direction, cal- 
culated from modelfit components within 7 mas from the 
co re, in our 2006 observat ions to the 2000 observations 
by Zavala & Taylor (2005), we see a change of 1 0°. We 
also do not dete cfas high positive RM values as |Zavala| 
fc Taylor| (|2005[), who see values up to 2000 rad m~^. 
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while our maximum values are near 500 rad m^^. The 
maximum gradient is detected about 3-7 mas from the 
core, where the RM changes from +500 rad m~^ to —600 
rad m~^. Further down the jet the gradient becomes less 
pronounced and we also detect less polarized jet emission. 



Asada et al. ( 2002| suggested that the gradient orig- 
inates from a helical magnetic field in the jet. Based 



on the large RM values observed in their study, [Zavala 
& Taylor ( 2005 ) preferred an external origin, possibly a 
sheath around the jet. The main argument they used was 
that internal Faraday rotation values of 2000 rad m~^ 
should cause severe depolarization in a uniform magnetic 
field, which was not observed. However, combinations of 
different magnetic field configuration and number of lines 
of sight can possibly explain high RM and complex po 



larization structure (see Sect. 4.2). Asada et al. (2008a I 



report variations in the transverse gradient ot '6C 273 be 
tween their observations in 1995 and 2002. They present 
several calculations ruling out the narrow line region as 
the origin of the Faraday rotation due to the variabil- 
ity and suggest the variations are caused by the external 
slower moving sheath changing over time scale of several 
years. Our observations are only three months apart, but 
there are still differences between the maps, especially in 
the region 2-5 mas from the core on the South side of 
the jet as can be seen in the values of slice 1. In the first 
epoch the RM changes from -1-450 to —300 rad m~^ but 
in our second epoch the change is from -1-500 to —760 
rad m^^. The large discrepancy in the negative values 
is also seen when examining component 14 in Table [3] 
located 2.7 mas from the core in Fig. |8] The component 
has not moved more than 0.2 mas (0.6 pc) over the two 
epochs, and still the values differ by over 400 rad m^^ 
so the change cannot be caused by the component illu- 
minating a different part of the Faraday screen. If the 
scale of RM variations in the screen were this small, we 
would not expect to see consistent RM values over the 
jet or well-defined gradients. 

This component also shows clear depolarization be- 
tween 15 and 8.1 GHz as the fractional polarization drops 
from 9.2 to 2.8% in the first epoch and from 8.4 to 1.4% in 
the second. This component is not on our list of isolated 
jet components due to the proximity of the brighter (in 
total intensity) component 12 and is therefore not shown 
in Fig. [T] Howev er, if we fit the polarization data sim- 
ilarly as in Sect, 
and 



4.2 



we obtain values — y |&| of —843 
-1001 for the hrst and second epoch, respectively. 
This would make the component the most depolarized 
in our sample. It is very difficult to explain such fast 
variations with external Faraday rotation, and therefore 
it is possible that we are seeing internal Faraday rotation 
in this case. Another alternative is that the variations 
we observe in the RM are due to interaction of the jet 
with a sheath. Chen ( 2005 ) observed variations of com- 
parable magnitude over similar time scales in 3C 273. 
He proposes that the fast variations could be caused by 
expansion of the components compressing the surround- 
ing medium increasing the magnetic field and electron 
density. However, it is difficult to explain the complex 
depolarization observations with this model. 

The observational signatures of large-scale helical mag- 
netic fields were recently studied from a theoretical per- 



spective by Clausen-Brown et al. (2011). They suggest 
that the best way to distinguish signatures of helical 
fields from interaction with external medium is to look 
for correlated behavior in the total intensity, spectral in- 
dex, and polarization profiles. In their model, the to- 
tal intensity, polarization and spectral index should have 
skewed profiles, so that a tail in total intensity is found 
on the same side of the jet where the polarization is lower 
and where the spectral index is steeper. The skewness of 
the profiles depends on the Lorentz factor and viewing 
angle of the jet. We have studied the total intensity, po- 
larization and spectral index profiles at the locations of 
the gradients in 3C 273 and show them for slice 2 in our 
March 2006 epoch in Fig. [9] It is ob vious that the profiles 
are sk ewed in the way predicted by Clausen-Brown et al 



(2011 1, supporting models with helical magnetic fields in 



3C 273. The spectral inde x gradient was also detected 
by Savolainen et al. ( 2008 1 . Unfortunately, such skewed 
signatures are m general difficult to detect due to beam 
effects and large errors in polarization towards the jet 
edges, and even in our observations the signature is not 
as clear in all jet locations. Our higher resolution and 
more sensitive VLBA follow-up observations will enable 
us to study this further. 

Another source which shows a significant transverse 
gradient in its March 2006 epoch is 3C 454.3, shown in 
Fig. [lOj The gradient is seen between 1-3 mas from 
the core and exceeds 3a. The magnitude of the gradient 
varies slightly depending on the chosen location with a 
max imu m of about 63 rad m~^mas~^. In the slice of 
Fig. 10 it is about 57 rad m~^mas~^ when the jet is 3 
beams wide. In our second epoch in June 2006, the gra- 
dient is not as clear, but that can be attributed to lower 
data quality in the 8.1 GHz band during that epoch, as 
a smaller region of the jet is visible above the noise level. 
Another complication arises from the bending of the jet 
because it is difficult to determine the transverse direc- 
tion when the jet bends. In 3C 454.3 we have chosen the 
local jet direction when studying the gradient, but it is 
obvious that the gradient is no longer seen further down 
in the jet after it bends. We see variations in the RMs 
of the jet components 1 and 2 which are 8.6 and 6.1 mas 
from the core, respectively, as well as inverse depolariza- 
tion in several components, and this could point towards 
internal Faraday rotation as seen i n 3C 273. Interest- 
ingly, Fig. 27 in [Zavala &: Taylor] ( |2003[ ) seems to hint 
at a RM gradient in the same direction, alt hough this 
was not reported by Zavala & Taylor (2003), who were 
not concerned with the possible presence of transverse 
gradients in their RM maps. When our observations are 
combined with total intensity, polarization and spectral 
index observations (Zamaninasab et al. 2012, in prep^ 
they s eem to follow a modification of the Clau sen-Brown] 



et al. ( 2011, ) model. Details of this modeling will be pre- 
sented m Zamaninasab et al. 

Additionally, several other sources show interesting 
transverse RM structures. In 2230-1-114 we detect a gra- 



dient of 144 rad m mas 



at 3cr level about 7 mas from 
Based 



the core where the jet is 1.9 beams wide (Fig. 11 ) 
on our criteria above, this can be considered as a signifi- 
cant gradient, but it is more difficult to tie it to any spe- 
cific model because the region over which the gradient is 
detected is small. Therefore we have included the source 
in follow-up VLBA observations which are designed to 
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Fig. 8. — Rotation measure maps of 3C 273 with high-RM pixels in the cores blanked for our March (left) and June (right) epochs. The 
black lines show slices transverse to the jet from which the gradients (bottom 3 panels) have been taken. The beam size along the slice is 
shown in each plot with a scale bar. A color version is available in the online edition of the Journal. 



give better sensitivity and resolution in hope of confirm- 
ing the gradient and modehng it in more detaiL The 
results of the follow-up observations will be presented in 
a separate paper. 

Similarly, 0923-1-392 shows a significant total gradient 
of over 624 rad m^^ where the polarized jet is 2.6 beams 



wide (Fig. 11 1. The gradient, however, is confined within 



depends on the high RM region at the Northern side of 
the jet. Interestingly, the high RM values are seen right 
where the jet is shown to bend o utside the line-of-sight of 
the observer ( Alberdi et al. 2000) and therefore could be a 



sign of the jet interacting with the intergalactic medium. 
The polarization structure we observe is also consistent 



with the observations of Alberdi et al. ( 2000 1 where the 
change in polarization over several years is shown to be 
consistent with a moving component interacting with a 
stationary feature where the jet bends. Another alterna- 
tive could be that the three dimensional geometry of the 
jet is complex and the North and South side of the jet 
probe different regions along the jet so that the other is 
further downstream and we are seeing effects of different 
optical thickness. With the present data it is difficult to 
distinguish between the alternative models. Therefore 
this source was also included in our follow-up VLBA ob- 
servations. 

We do see hints of transverse gradients in three other 
sources as well, but none of these fulfill all of our criteria. 
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Fig. 9.— Transverse slices of 3C 273 along slice 2 in Fig. |8] left 
panel for RM (top panel), total intensity (2nd panel from the top) 
linearly polarized intensity (3rd panel from the top), and spectral 
index between 15.3 and 8.1GHz (bottom panel). 15.3 GHz obser- 
vations are shown in (red) filled circles and 8.1 GHz observations in 
(black) crosses. The rms error in total intensity is 1.2 mjy/beam 
and 4.0 mJy/beam at 15.3 and 8.1GHz, respectively. The rms 
error in linear polarization is 0.6 and 0.8 mJy/beam at 15.3 and 
8.1 GHz, respectively. The error in spectral index varies from 0.004 
in the center of the jet to 0.3 at the edges. The beam size along 
the slice is shown in each plot with a scale bar. Color version is 
available in the online edition of the Journal. 



For example, in 2134+004 (Fig.|Tjl39) the change in the 
RM is more than 2a, but it is detected only at a single 
location of the jet and is only 65 rad m~^ mas~^ where 
the jet is 1.8 beams wide, so we do not call it significant. 
In 0945+408 the RM map (Fig. [l]67) shows clearly two 
different RM regions, and even though the jet is more 
than 2 beams wide, the gradient is within 2a errors. A 
similar gradient is seen in 1641+399 (Fig. [1]110) where 
it also is within 2a errors. Therefore we cannot call the 
gradients in these sources significant but have included 
them in the follow-up higher sensitivity observations for 
further study. 

We have also compared our maps to other studies re- 
porting transverse RM grad ients in the sources in our 
sampl e. Asada et al. (2008b) and Reichstein & Gabuzda 
(2011) report a transverse gradient of a tew hundred 
rad m~^ in the source 0333+321 in VLBA observations 
between 5 and 8 GHz. We do see a gradient of simi 



lar magnitude (Fig. [T]21), but it extends only one beam 
width across the jet, and it is also very much dependent 
on low S/N-jet edge pixels. Additionally, it does not ex- 



tend over the whole length of the jet. Some of this may 
be due to our higher observing frequency which causes us 
to detect less polarized emission across the jet. Therefore 
we do not consider this a robust gradient in our images 
but merely suggestive of a possible gradient. IReichstem] 
fc Gabuzda, (2011, ) also report gradients in both the core 
and jet of 1150+812 but in both cases the slices they 
take are less than two beams wide at the location of the 
gradient. We do see a change in RM values in a similar 
direction (Fig. [l] 79) but again it is very much dependent 
on unreliable edge pixels and also the jet is only one beam 
width across at our 8 GHz resolution so that we do not 
consider the gradient to be robust. 

In 3C 120 (Fig. [l}28) we do not detect the transverse 
RM gradient that was seen by Gomez et al. (2008) in 
observations made in 20 01 at 15, 22 and 43 GHz. I'his it 
not unexpected because Gomez et al. (2011 ) demonstrate 
in their Fig. 9 how a different region oi the jet in 3C 120 is 
seen in their 2001 and 2007 observations at 15 GHz. Our 
observations are close to their latter epoch, and therefore 
we also see a different part of the jet and do not detect 
the gradient. 

Other repo r ted so urces include 0836+710, for which 
Asada et al. (2010) find a possible gradient of ~150 
rad m~"^ in observations between 4.6 and 8.6 GHz. The 
gradient is at a location in which we were not able to ob- 
tain good A^-fits (Fig.[l]59) and therefore we cannot con- 
firm the result. In any case, the jet is not well-resolved 
transversely in our observa tions. Another source of inter- 



est is 1803+784, for which [Mahmud et al.| ( |2009[ ) report 
RM gradients in four epochs observed with several differ- 
ent frequency setups from 4.6 to 43 GHz. In our map we 
do not see any indication of a gradient (Fig.[l|l22) , but it 
must be noted that our map does not extena as far down 
the jet as some of theirs. This is because in their lower 
frequency maps they detect polarization fu rther down 
the jet compared to our 15 GHz observations. lO'Sullivan] 



& Gabuzda (20091 studied the Faraday rotation in six 
blazars observed m July 2006 between 4.6 and 43 GHz, 
and detected a transverse gradient in two of them in mul- 
tiple frequency ranges. In 0954+658 they detect a gra- 
dient which is mainly visible in the low-frequency RM 
map from 4.6 to 15.4 GHz. Even in this case the jet is 
unresolved. In the 7.9 - 15.4 GHz RM map, very close 
to our resolution, they see an indication of the gradient, 
but again the jet is unresolved. Our map (Fig. [l]69) 
shows similar changes in the RM from positive to near 
rad m~^ at the same location, but the jet is less than two 
beam widths wide and the gradient is also very much de- 
pendent on the edge pixels, which we have noted in the 
previous section to be unreliable. We cannot compare 
the core RM values because we do not find good A^-fits 
in the core region. We also detect changes in the jet RM 
because in our observations the jet RM varies between 
an d 600 rad m~^ while in t he sam e frequency range map 
in |0'Sullivan & Gabuzda] (|2009| the values are between 

T 



and 200 rad m~^. This means that we are either seeing 
a different part of the jet, or the RM has changed over a 
time span of three months. The gradient they detect in 
1156+295 is highly dependent on a few unrehable edge 
pixels. Additionally it is within one beam width from 
the core a nd therefore unreliable as shown by the sim- 
ulations of Broderick & McKinney (2010). We do not 



Faraday rotation in parsec-scale AGN jets 



17 




2251+158 
3006_06_15 




-5 -10 -15 

Relative E.A. (mas) 



-5 -10 
Relative E.A. (mas) 





distance along the slice (mas) 



distance along the slice (mas) 



Fig. 10. — Rotation measure maps of 3C 454.3 with high-RM pixels blanked for our March (left) and June (right) epochs. The beam 
size along the slice is shown in each plot with a scale bar. Color version is available in the online edition of the Journal. 



detect any polarized jet emission in this source at our 
epoch (Fig.jl]79) and therefore do not search for gradi- 
ents. We do not obtain good A^-fits in the core but the 
few pixels where we obtain acceptable fits agree with the 
R M values of | ; Sullivan & Gabuzda (2009). 

'Contop oulos~et af.| ([2009) list 29 sources for which a 
gradient can be seen in RM maps collected from the lit- 
erature. Our sample includes 25 of these a nd we have 
looked for gr adients in the 18 sources in which [C"ontopou-] 
los et al. (2009) report a gradient at least 2 mas from the 
core. I'hese include 3C 273, 3C 454.3 and 2230-^114 and 
two others which are in our list of sources showing hints 
of gradients. In the remaining 14 we either do not detect 
RM outside a beam width from the core (5 sources) or 
do not see any indication of a gradient in the jet at our 
epochs. We note that some of these observations were 
done at longer wavelengths and therefore may have been 
more sensitive to polarization and Faraday rotation. 

Based on our results of a very large sample of RM maps 
we conclude that robust mas-scale RM gradients are dif- 
ficult to detect at > 3cr level in the inner jets of AGN. 
The major cause for this is that the majority of the jets 
we observe are not transversely well-resolved and there- 
fore we do not consider the gradients as ro bust based on 
our detailed simulations in Appendix |C.2| We note that 
multiepoch observations may help to confirm other gra- 
dients because in principle, detections of three 2cr gra- 



dients in the same source at the same location would 
give an overall significance of 3cr. However, this requires 
that the significance of the individual detections is deter- 
mined using statistically correct error bars derived from 
error propagation of Q and U rms error and accounting 
for additional D-term and CLEAN errors as described in 
AppendixlBl In addition to the four sources where we de- 
tect significant transverse gradients, our sample includes 
only five other sources which have polarized jets wider 
than two beam widths and show no sign of a transverse 
gradient. The reduced sensitivity and angular resolution 
makes the detection of gradients very difficult in objects 
at higher redshifts. For example, if we restore the maps 
of 3C 273 with a beam size corresponding to angular res- 
olution of a z = 0.5 object and reduce the flux density by 
a factor of 10 to achieve a typical flux density of a source 
at z = 0.5, the gradient would no longer be significant. 
Higher sensitivity and better angular resolution observa- 
tions would help to solve the problem. Given the existing 
constraints of the current VLBI arrays, this is a strong 
science motivation for high-sensitivity space- VLBI. 

4.4. Time variability 

As was noted in the previous section, there is clearly 
variability in the RM values over time scales of years or 
even months. If the Faraday rotation is caused by an 
external screen very far away from the jet, for example 
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Fig. 11. — RM maps and transverse RM slices of 2230+114 (left) and 0923+392 (right). Beam size along the slice is shown in each plot 
with a scale bar. Color version is available in the online edition of the Journal. 



in intergalactic clouds or even in the narrow line region 
of the AGN, the RMs observed in fixed locations (with 
respect to the core) in the jets of AGN should remain 
constant over time scales of years. In order for external 
clouds to cause variations over time scales of years, they 
should be so small in size that the number density is 
not large enough to cause significant variations in the 
observed RM. Therefore it is important to verify that 
the variations are seen in the same parts of the jet and 



are not due to changes in the jet position angle ( Gomez 



Twenty sources in our sample were observed twice dur- 
ing the 12 month period, and in ten AGNs we measured 
RM values in both of the epochs. Additionally, we can 
co mpare our observat ions to the 8-15 GHz observations 
of |Taylor| ( pOSj j2000t and [Zavala fc Taylor] p002l [2003 
' 20U4p wliere RM maps for 36 sources m our sample, 
mamly quasars, are shown (3C 279 is shown in three of 
them). These observations were obtained between 1997 
and 2000, and our comparisons are therefore affected by 
the long gap between our observations. Another factor 
that will affect at least some sources is that Zavala & 
Taylor did not remove the contribution of Galactic Fara- 
day rotation from the ir maps . This affects, for exam- 



ple, 2200+420 where [Zavala fc Taylor[ ( |2003[ ) observe 
a jet RM of —287 rad m~^ which will reduce to —136 
rad m~^ if the Galactic Faraday rotation is accounted 



for. This value is consistent within our error bars with 
our jet RM values of about rad m~^. RM maps of five 
of our sources wer e obtained in July 2006 by [0'Sullivan| 



fc Gabuzda ( 2009[) which allows us to probe shorter time 
scale variations. I'hey include RM maps between 7.9 and 
15 GHz for most of their sources, close to our resolution. 
Out of these five, 0954+658 and 1156+295 were already 
discussed in the previous section. 

4.4.1. Variability in the core RMs 

In four sources which were observed twice during our 
program (0215-f015, 0716-1-714, 0834-210, 0847-120), 
we detect RM values only in small regions or near the 
core and the values are all within la error bars in the 
two epochs. In 0219+428 and 2200+420 our observa- 
tions agree very well. We do not detect good A^-fits in 
the cores of these two sources in either of the epochs. In 
th e case of 2200+420 this agrees well with observations 
of [Mutel etal] ( |2005[) who observed 2200-^420 over nine 
epochs between lU^ and 2002 at 15, 22 and 43 GHz, 
and did not detect good A^-fits in four of their epochs. 
They showed that this is due to blending of multiple com- 
ponents in the core region. The components were seen 
as separated in their 22 and 43 GHz observations but 
were blending together when restored with the 15 GHz 
beam. This could also be the case in 1418-1-546, where 
four Gaussian components are required to fit the region 
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2 mas from the core at 15 GHz. lO'Sullivan fc Gabuzdal 
( |2009 ) also observed this source in July 2006 between 
4.6 and 43 GHz. In our February 2006 epoch we obtain a 
similar core RM value as they do, although we do not de- 
tect the positive RM patch as seen in their 7.9-15.3 GHz 
map. In our November 2006 map the core values change 
by 200 rad m~^, compared to our previous epoch and 
their map from July. In this epoch, we also detect the 
p ositive RM patch seen in their map. 

O'Sullivan fc Gabuzda] ([2009^ observed 1749+096 in 



July 2006 but they do not obtain any good fits, which 
is surprising because we detect RM values in the core of 
1749+096 in our June epoch, observed only two weeks 
earlier. They attribute the inability to obtain good fits 
to a possible flare in the source which may affect the 
polarization structure. Indeed, their peak flux in total 
intensity in the 7.9 GHz band is 4.2 Jy/beam while in 
our 8.1 GHz observation it is 3.5 Jy/beam so that a flare 
is probably on-going and could affect their July observa- 
tions. In fact, a flare peaking in October 2006 is seen in 
the 15 GHz MOJAVE observation^ Additionally, our 
wavelength coverage is not as wide as theirs so it is pos- 
sible that we do not detect the non-A^ behavior they see. 
Variations in the core RMs ar e also seen in compa rison 



2000) and 



Only in 2230+114 do we see differing values in the jet 
RM, even though it looks like we a re looking at the same 
portions of the jet. |Taylor| ( |2000[ ) detect a RM of -185 
rad in the jet of 2230+114 (component C in their 
paper) which is almost at the same location as our com- 
ponent 3 with RM of +173 rad m~^. The different sign 
in the observed RMs indicates a change of direction in 
the line-of-sight component of the magnetic field. It will 
be interesting to compare these to our follow-up obser- 
vations to determine if the sign has remained the same. 

Even though at first glance it looks like the RM val- 
ues in the jets of AGN change over time scales of several 
years, a detailed comparison taking the differing loca- 
tions of the polarized components into account shows 
that we do not detect significant variations in the jet 
RMs of the majority of AGN. The fact that we detect 
non-zero RMs in the jets even after the Galactic Faraday 
rotation contribution is taken into account suggests that 
the RM is occurring outside our own Galaxy, for exam- 
ple in the na rrow line region of the AGN (e.g. Zavala & 
Taylor|2003| ). In three sources which show signs of trans 



of our observatio i is to those o f Taylor ( 1998 
Zavala & Taylor (2003 2004). In general, the core val- 
ues between our maps and theirs differ significantly and 
there are only five sources where we detect similar core 
RM values. This can also be attributed to the blending 
of components within the finite beam. Alternatively, if 
the Faraday rotation is internal to the jet, changes in the 
particle density or magnetic field strength due to newly 
emerging components could cause the variations. In six 
sources the core RM in our maps has a different sign 
than in their observations which indicates a change in 
the direction of the line-of-sight component of the mag- 
netic field. However, due to the long gap between our 
observations and the complexity of the core regions in 
these sources, it is difficult to distinguish if this is due 
to a sign-reversal in an ordered magnetic field or due to 
changes in properties of new components. 

4.4.2. Variability in the jet RMs 

3C 273 and 3C 454.3 were seen to have large variations 
in the jet RMs over time scales of three month in our 
observations, which was discussed in the previous section. 
In some sources it is clear that we are seeing a different 
part of the jet at successive epochs and sometimes we do 
not detect the polarized jet emission at all, or vice versa. 
This is seen, for example, in complex sources such as 
1253—055 and 1418+456 where we do not see the exact 
same part of the jet because the same jet components are 
not detected or the components have moved significantly 
between the two epochs. 

In many sources we still see fairly similar jet RM val- 
ues within the error bars. In the jet of 2200+420 we 
observe R M values con sistent with rad m~^, in agree- 
ment with Mutcl ct al. ( 2005 ) . This source was observed 
by [O^Sullivan fc Gabuzda^ (2009) in July 2006 between 
4.6 and 43 GHz where they obtain slightly higher RM 
values in the jet than we do. However, these are still for 
the most part within the error bars of our observations. 

^ http:/ /www. physics. purdue.edu/astro/MOJAVE/sourcepages/1749 



verse rotation measure gradients, 3C 273, 2230+114 and 
3C 454.3, RM time variations are seen in the jet as well. 
In 3C 273 and 3C 454.3 the variations happen over time 
scales of three months, pointing to either internal Fara- 
day rotation or interaction between the jet and the Fara- 
day rotating material, possibly a s heat h around the jet, 
as discussed in more detail in Sect. 14.31 If the rotation is 
internal, the rotating plasma is either the low-energy end 
of the relativistic electron population or thermal plasma 
intermixed with the emitting plasma, and therefore fast 
changes in the magnetic field of the jet can cause fast 
variations in the observed RM. If the jet is surrounded by 
a sheath, the sheath must be mildly relativistic because 
we do not detect counter-jet emission in these sources. 
Therefore over time scales of years, it is possible that 
the RM v ariations are caused by changes in the screen 
itself (e.g. Asada et al. 2008a). This could explain the 
variability in 2230+114 but not in the two other sources 
where the gas in the sheath would not have moved suffi- 
ciently over the time scale of three months to cause sig- 
nificant variations. In sources with fast apparent super- 
luminal motion, it may be possible to observe variations 
over time scales of months when the components illumi- 
nate a different part of the Faraday screen (e.g. Gomez 
|et al.||201l] ). In these sources the same RM should be 
seen when another component passes the exactly same 
part of the jet which could then be studied with frequent 
multiepoch observations. This seems not to be the case 
in 3C 273 or 3C 454.3 where the components have only 
moved about 0.1 mas between the two epochs. We note 
that if the variations happen only in small parts of the 
jets, our ability to detect variability is easier for nearby 
sources with resolved jets. 

4.5. Comparison to simulations 



Re cently , usin g GRMHD simulations, Broderick fc 

|McKinney ( 2010 ) produced simulated RM maps for AGN 
jets with large scale toroidally dominated magnetic fields. 
They produced maps with different resolutions and beam 
sizes to additionally study the effect of finite beam size 
on RM observations. Their results showed that any RM 
"^|MH5fetft?i observed within a beam width from the core 
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Fig. 12. — Maps of the rotation measure sign in nine example objects showing bilateral structure (see text for details). Dark gray (blue) 
areas are positive RM and light gray (red) areas negative RM. Color version is available in the online edition of the Journal. 



are unreliable at resolutions below 43 GHz at VLBA. Fur- 
ther down the jet it should be possible to detect true RM 
gradients if the magnetic field is toroidal. Formally, with 
infinite sensitivity, even in unresolved jets it may be pos- 
sible to detect gradients of the right sign, with the mag- 
nitude severely suppressed. They did not, however, take 
the effect of noise into accoun t in t heir simulations, which 
we have shown in Appendix |C.2| to have a major effect 
on the reliability of gradients in unresolved jets. They 
also showed that if the toroidal component of the field 
is made to vanish, any significant gradients will vanish 
too. They conclude that linear, resolved transverse gra- 
dients are due to large-scale toroidally dominated mag- 
netic fields within the Faraday screen. This provides us 
a way to compare our observations with their simula- 
tions. It is diff icult to make quan titative comparisons, 
but IBroderick & McKinney (20101 show that even maps 
of the KM sign should show the bilateral structure if 
the toroidal magnetic field component dominates. Their 
Fig. 12 shows examples of sign maps for different res- 
olutions and other parameters, such as the jet viewing 
angle and black hole mass. Additionally they present a 



model in which the toroidal field is made to vanish, and 
a random foreground Faraday screen is inserted. 

In nine sources visual inspection of the sign maps re- 
veals clear bilateral structures as shown in Fig. [12]. These 
include 3C 273, 3C 454.3, 0923-1-392 and 2230-fll4 in 
which the gradie nt is significant based on our detailed 
analysis in Sect. 4.3 In two other sources (0945-1-408 
and 1641-1-399), in we which detect the gradient at a 2<j 
level, the sign maps show different signs on both sides 
of the jet along the whole polarized region of the jet. 
These are candidate sources for toroidally dominated 
large-scale magnetic fields. However, one should note 
that for example in the case of 1611-1-343, even though 
the sign map shows a bilateral structure, any possible 
gradient is within the Icr error bars, and the jet is over 
four beams wide at its widest location. In 0333-1-321 and 
0754-1-100 the jets are not sufficiently resolved for us to 
make any conclusions about the significance of the gradi- 
ent. If the jets are not viewed directly from the side in the 
jet frame, it may be that we are seeing only the positive 
or negative side of the possible toroidal field. Therefore 
we have produced sign maps where the median jet region 
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Fig. 13. — Maps of the rotation measure sign in nine example objects that represent the typical sources in our sample (see text for 
details). Dark gray (blue) areas are positive RM and light gray (red) areas negative RM. Color version is available in the online edition of 
the Journal. 



RM is subtracted from the map to see if we can detect 
more sources with bilateral structures. The re are only 



two such sources, 1150+812 and 2005+403. Reichstein 



& Gabuzda (20121 report a sign change centered around 
zero in 1150+812 in lower frequency data between 4.6 
and 8.9 GHz, while in our map the center is about 240 
rad m^^. In our map the change in magnitude also hap- 
pens nearly transverse to the jet while in their map it is 
more along the jet. The jet of 1150+812 is only one beam 
width wide in our observations and therefore we cannot 
make conclusions about the significance of the gradient. 
In 2005+403 the jet is 1.5 beams wide in polarization and 
the values are centered around ~110 rad m~^, but are 
all within one a of each other and therefore the gradient 
is not sig nific ant. 

In Fig. [13] wc show 9 example cases that we believe 
represent typical types of behavior in the sign maps. In 
many sources we see how the RM sign changes along the 
jet direction (e.g., 0003-066 and 1222+216). These kind 
of structures look similar to simulations of unresolved jets 
with zero toroidal magnetic field component or random 
magnetic field. Multi-epoch studies of changes in the RM 



and polarization structure could give more insights into 
the nature of the screen and if it is internal or external to 
the jet. In the bottom row we show examples of sources 
in which the jet is transversely resolved but the sign maps 
show a very complex structure similar to the simulated 
case with resolved jet and zero toroidal magnetic field or 
random foreground screens (e.g., 1253—055) or no change 
in the sign of the RM. The majority of t he sources in our 
sample show sign maps similar to Fig. [13] or sign maps 
without any change in the sign. We note, however, that 
the majority of our jets are transversely unresolved and 
therefore it may be possible that we simply cannot detect 
the bilateral structure in our maps. 

5. CONCLUSIONS 

We have observed a sample of 191 AGN jets (20 at two 
epochs) at 8.1, 8.4, 12.1, and 15 GHz with the VLBA to 
determine the Faraday rotation measures in the parsec- 
scale jets of these sources. One motivation for the study 
was to find out how much Faraday rotation would af- 
fect the EVPAs of the MOJAVE program observations at 
15 GHz. The polarization is sufficiently strong to study 
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Faraday rotation in 149 sources (159 maps), and in over 
80% of the sources the median observed RM is less than 
400 rad m~^, which would rotate the 15 GHz EVPAs by 
less than 10 degrees. Additionally, we do not see signif- 
icant variations in the jet RMs of most of the sources 
so that these results can be extrapolated to other MO- 
JAVE epochs as well. However, there are several sources 
with high observed RMs of larger than 1000 rad m~^ 
so that when studying polarization in individual sources 
(especially at frequencies lower than 15 GHz), one needs 
to take the Faraday rotation into account. The highest 
intrinsic redshift-corrected absolute RM 8975 rad m~^ is 
detected in a jet component 0.5 mas from the core in the 
high-redshift quasar 0642+449 at ^ = 3.396. 

We find the quasars to have generally higher abso- 
lute RMs than the BL Lac objects (median |RMobs| 144 
rad m~^ vs. 79 rad m~^) and the core components in 
the sources to have higher RMs (median |RMobs| 171 
rad m^^) compared to the jets (median |RMobs| 125 
rad m~^) which is also seen in both quasars and BL Lacs 
separately. In quasars, we also detect a significant neg- 
ative correlation with the magnitude of the RM and the 
de-projected distance of the component from the core. 

We perform detailed simulations of the measurement 
error of RM, and focus especially on the effect of noise 
and finite restoring beam size on putative transverse RM 
gradients. Because of the finite beam size, neighboring 
pixels are not independent and this has to be taken into 
account when interpreting radio interferometric images. 
For a typical jet there are only few independent measure- 
ments across a jet of apparently many pixels. Our simu- 
lations show that it is possible to obtain spurious, yet sig- 
nificant looking transverse RM gradients if the polarized 
region in the jet is less than two beams wide. We give 
several guidelines on how the errors in RM maps need 
to be taken into account to determine the reliability of 
a transverse gradient. The reliability is strongly depen- 
dent on the jet width, and the jet should preferably be 
two beams wide and more than 3cr in significance to call 
the gradient reliable. The errors in the RM maps should 
be calculated from the variance-covariance matrix of the 
least-squares fit where the EVPA errors are determined 
using error propagation from Stokes Q and U rms errors, 
additionally accounting for effects due to the CLEAN 
procedure. Following these guidelines, we detect signifi- 
cant transverse RM gradients in the flat spectrum radio 
quasars 3C 273, 3C 454.3, 0923-^392 and 2230-fll4. In 
3C 273 the RM is for the first time seen to change sign 
along the transverse slice, giving further support for a 
helical magnetic field in the jet. The main reason why 
we are not able detect gradients in more sources is that 
the jets are insufficiently resolved at our lowest observing 
frequency of 8.1GHz. In addition to these four sources, 
there are only five which have wide-enough polarized jets 
where a gradient is not detected. The reduced sensitivity 
and angular resolution makes the detection of gradients 
very difficult in objects at higher redshifts. Higher an- 
gular resolution observations of much greater sensitivity 
are needed to study this phenomenon further. 

Comparison of our RM maps to earlier studies con- 
firm that significant temporal RM variations in the core 
components of AGNs are common. This could due to 
multiple polarized components blending within the finite 



beam at different times or intrinsic changes in the mag- 
netic field or particle density if the rotation is internal 
to the jet. Higher resolution observations at higher fre- 
quencies are required to uncover the true RMs in the 
cores of many of the sources. In almost all the cases 
where our jet RM values differ from earlier studies, it 
can be explained with different part of the jet being illu- 
minated by the components at different times. Only in 
3C 273, 3C 454.3, and 2230+114 do we detect variations 
which cannot be explained by moving components. In 
2230+114 the variations happen on time scales of years, 
which could be possible if the Faraday rotation occurs 
in a mildly-relativistic sheath around the jet. In 3C 273 
and 3C 454.3 we see variations over a time scale of three 
months which are difficult to explain with external Fara- 
day rotation models, so more likely the rotation is inter- 
nal to the jet. This is further supported by observations 
of inverse depolarization structures in the jet components 
of these two sources which show the 8 GHz polarization 
to be higher than at 15 GHz. This is the opposite trend 
with wavelength expected from standard external depo- 
larization and may be explained by combining internal 
Faraday rotatio n with helical or loosely tangled random 
magnetic fields ( Homan|2012 ). However, in the majority 
of isolated jet components we studied, the relationship 
between depolarization and RM can be explained with a 
purely random external Faraday screen viewed through 
a small number of lines of sight. In these sources the 
rotation could be caused by intergalactic clouds or the 
narrow line region of the AGN. 
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APPENDIX 

A. EFFECT OF IMAGE ALIGNMENT ON RM MAPS 

To study the effect of image alignment errors on the 
rotation measure maps, we used 12 sources for which 
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we have a very good correspondence between the shifts 
from our 2D cross correlation and from opticahy thin 
component positions. In all these sources the difference 
between the two shifts is less than 0.1 mas. We assumed 
the shifts in these cases to be correct and introduced ar- 
tificial shifts between the maps to study the effect on 
both spectral index and rotation measure maps. The 
shifts were applied to the maps of all the other frequency 
bands with respect to the 15.3 GHz maps. In the case 
of spectral index maps, the effect was large, with some 
sources showing clear differences with shifts as small as 
0.03 mas. In other sources, even a shift of 0.18 mas, al- 
most twice our median shift for the whole sample, did 
not cause major changes in the spectral index maps. For 
each of the test sources we calculated rotation measure 
maps for shifts that were seen to create spurious features 
in the spectral index maps. The effect on the rotation 
measure maps was much smaller, and shifts as large as 
0.18 mas did not affect the general structure. In most 
of the cases small differences could be seen in the edge 
pixels or close to patches of low signal to noise, which in 
any case are considered unreliable. In only two maps did 
we see appreciable differences in the jet rotation mea- 
sure, and in these cases the wrong shifts were 0.09 and 
0.15 mas, but those were also in very complicated regions 
of the jet. Even in these two cases the general structure of 
the RM map did not change and our conclusions would 
be unaffected. Therefore we conclude that even if our 
image alignment is off by 0.15 mas, it should not affect 
the results from our rotation measure maps, especially 
as we are not using the edge or low signal-to-noise re- 
gions to make conclusions about the RM structure. By 
using the spectral index map as an additional indicator 
of the goodness of the alignment we ensure that our ro- 
tation measure maps arc not affected by false alignment 
between the different bands. 

An example of spectral index maps in the ca se w here 
the alignment is very important is shown in Fig. 14 The 
figure shows spectral index maps without and witn shift- 
ing between all the bands. The correct shifts were ob- 
tained using the 2D-cross correlation and are 0.39 and 
0.03 mas in right ascension and declination between the 
15 and 12 GHz bands, 0.45 and 0.03 mas between the 
15 and 8.4 GHz bands and 0.33 and 0.03 mas between 
the 15 and 8.1 GHz bands, and are shown in the bottom 
panel of Fig. [T4j Analogous to results of 'Kovalev et al.' 
(2008), it is clear that the maps without shifting that 
have b een aligned based on the phase center (top panel 



of Fig. 14) have spurious features, such as optically thin 



regions apparently upstream of the core and a spectral 
index gradient transverse to the jet, which all disappear 
when the images are properly aligned. 



B. 



ERRORS IN POLARIZED FLUX DENSITY, 
FRACTIONAL POLARIZATION AND EVPA 



Errors in polarized flux density Up are typically as- 
sumed to be the average rms error in the Q and U im- 
ages. Additional sources of error are errors due to the 
CLE AN procedure whi ch can be difficult to overcome 



(e.g.. Lister et al. 2001) and errors due to instrumental 



polarization (U-te rms) which are not evenly distrib uted 



across the images ( [Roberts, Wardle fc Brown||1994[ ). We 
used simulations to examine whether the error estimates 



we use are consistent with the above error contributions 
to the data. The simulations were carried out in several 
steps: 

1. A Stokes I model of the source was created using 
calibrated {u,v) data of a real source and CLEAN 
in Difmap. 

2. Stokes Q and U models were created by setting Q 
and U to be known fractions of Stokes I for each 
CLEAN component. Without any noise added, 
this corresponds to a uniform fractional polariza- 
tion and EVPA across the source. 

3. The original {u,v) data were loaded into AIPS and 
the task UVMOD was used to replace the real data 
with the values produced in the previous step. Ad- 
ditionally, random noise of the same order as seen 
in our real data was added. 

4. The UVMOD task was repeated 100 times to pro- 
duce 100 simulated {u,v) data sets with random 
noise added. 



5. 



In the case of D-term simulations, additional ran- 
dom D-term error with a standard deviation of 
0.002, determined from the scatter in the D-terms 



of our data (see Section 2.2 for details on the de- 
termination of the stable D-tcrms), was added to 
each set using the task SPLIT. 

6. The simulated {u,v) data were then imaged in 
Difmap following the same procedure as for the real 
data to obtain 100 images in Stokes I, Q and U. 

7. The rms in each image was obtained by shifting the 
map by 1 arcsec and calculating the rms using the 
command 'imstat' in Difmap. 

The process was repeated for two sources, 0333+321 
and 1226-1-023, to verify that the errors are general and 
do not depend on the specific structure of a given source. 
In order to address the errors due to the CLEAN pro- 
cedure and the clipp ing of the {u,v) coverage of our real 
data (see Sect. 2.1), the simulations were repeated for 
the original un-clipped 15 GHz data and the data at all 
the frequency bands were treated with the same cutoffs 
as our real images. We then studied the distributions of 
polarized flux density, fractional polariza tion an d EV PA 
in individual jet locations shown in Figs. [15] and [THj We 
chose locations at both the brightest parts of the jet and 
at the jet edges to see how well the error formula repro- 
duces the standard deviation of the simulations. 

For each location we determine the expected value and 
its expected error for all the parameters. For I, Q and U 
images the expected values are determined from low noise 
images which have been created using the same simula- 
tion procedure but adding only 1% of the typical noise in 
the task UVMOD. Additionally, the last CLEAN step (6) 
was replaced by a simple restoration of the model com- 
ponents from step 1 to have comparison images which do 
not suffer from CLEAN errors. The expected rms value 
is taken from the simulated maps with noise added at a 
location 1 arcsec from the map center. The expected p, 
EVPA and m values are calculated from the low noise 
images with p = ^/Q^TTP, EVPA = 1/2 tan" \U/Q) 
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Fig. 14. — Spectral index maps between 15 and 12 GHz (a and d), 15 and 8.4 GHz (b and e), 15 and 8.1 GHz (c and f) aligned by map 
center (top panel a-c) and by using correct shifts from the 2D-cross correlation (bottom panel d-f). Contours correspond to the total 
intensity at 15 GHz. A color version is available in the online edition of the Journal. 
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Fig. 15. — Locations of the simulated distributions in the jet of 
1226+023 overlaid on the 15 GHz total intensity contours. Star 
shows the component discussed in the error analysis. 



Fig. 16. — Locations of the simulated distributions in the jet of 
0333+321 overlaid on the 15 GHz total intensity contours. Star 
shows the component discussed in the error analysis. 



and m = p/I. The expected error values are calculated 
using the following equations: 



ctevpa 



and 



2(Q2 + [/2) 2p' 



(Bl) 
(B2) 
(B3) 



where aq and aij are the Q and U rms values in most 
of the simulations, and Q and U are the Stokes param- 
eters in the given pixel. In |B3| we assume that the ct/ 
term can be neglected as {p/ I'^f will be very small for 
all our components. For the simulations which take the 
instrumental polarization into account, an additional er- 
ror (TDterm defined as 

0-002 [-, 

V^ant * Nw * A'^^can ^ 

(B4) 
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where iVant is the number of antennas, Nip is the num- 
ber of IFs, A'scan is the number of scans with independent 
parahactic angles and /peak is the peak total intensity of 
the map ( [Roberts, Wardle fc Brown] [T994[ ). This ad- 
ditional error term is added m quadrature to the rms 
errors. The factor 0.002 in the equation is determined 
from the scatter of the D-terms in our data. In our data 
the number of antennas is 10, the number of IFs is 4 
for 15.4 and 12.1 GHz, 2 for 8.4 and 8.1GHz, and the 
number of independe nt scans is 4. The equation is de - 



fined and explained by Roberts, Wardle & Brown (1994), 
where they study the ettects ot mstrumental polarization 
with detailed simulations. For example, they show that 
the D- term errors scatter across the image (factor 0.3 
in Eq. B4) which seems to be supported by our simula- 
tions. I'herefore we include a contribution fro m th e total 
intensity peak and the current location in Eq. |B4[ 

The D-term error is strongly dependent on the Stokes 
I and is not distributed evenly across the images. There- 
fore its effect is largest on the bright locations of the 
source, especially near the core and it is not accounted for 
in the rms errors. To demonstrate the effect, in Fig. [TT] 
we show distributions of total intensity I, Stokes param- 
eters Q and U, polarized flux density p, EVPA and frac- 
tional polarization m for 1226-1-023 at point (shown 
by a star) of Fig. 15 In this simulation Q and U were 
set to be 0.0707 xT in each component. Without any 
noise added, this results in uniform fractional polariza- 
tion TO = 0.1 and EVPA = —22.5 degrees. Above each 
histogram we give the expected value for each parameter 
and its e xpected error value calculated using Eqs. |Bl[|B2[ 
and |B3| The top six panels show 100 simulations with- 
out adding the D-term error in step 5 in the simulations. 
As can be seen from the expected and observed values, 
the observed standard deviations of Q and U are slightly 
larger than the measured Qrms a-nd C/rms and the distri- 
butions are not peaking at the expected values, which 
can be attributed to the CLEAN errors. In the bottom 
six panels, the simulation is repeated with additional D- 
term noise added, whi ch is also included in the expected 
error value using Eq. |B4| Now the expected and ob- 
served rms values correspond very well to each other. In 
this bright location of the source, the contribution from 
D-term errors to the error in polarized flux density is 
twice as large as the rms error (1.1 mJy compared to 
0.56 mJy). The distributions of other components are 
not shown here but they follow the same pattern, except 
that the D-term error contribution diminishes with de- 
creasing Stokes I so that in most of the components of 
3C 273 the rms errors and D-term errors are of the same 
order. In fainter sources, such as 0333-1-321 the D-term 
errors are typically negligible in the jet components and 
smaller than rms errors even in the core . This effect is 
correctly accounted for by using Eq. |B4[ 

In order to study the effects of (M,?;)-clipping and the 
CLEAN procedure, we simulated the original 15 GHz 
data of 0333-1-321 and the (M,?;)-clipped data restored 
to the 8.1GHz beam. In these simulations we set Q = 
-0.031623 X I and U = -0.094869 x I, which resuhs 
in uniform fractional polarization to = 0.1 and EVPA 
—54.217 degrees. The expected values were determined 
as in the case of 1226-1-023. In Fig. [T8| top six panels we 
show the distributions for the original images and in the 
bottom six panels the distributions for the images which 



have a reduced (M,w)-coverage corresponding to the real 
data and restored with the 8.1GHz beam size. These 
are taken at location 6 (shown by a star) in Fig. 16 to 
show the effects on a faint edge region. Our equations for 
the error values work very well showing that the CLEAN 
procedure does not affect the standard deviation of the 
distributions, but there is a large offset in the peaks of 
the distributions in both simulations. This is also seen 
in the lower frequency simulations, for which we show an 
example at 12.1 and 8.1 GHz (8.4 GHz is nearly identical 
to 8.1GHz and therefore not shown) in Fig. \T9\ This 
is due to the reduced (M,w)-coverage, especially at 15.4 
and 12.1 GHz, and the CLEAN procedure. By looking 
at the offsets in all the component locations in the dif- 
ferent sources and all the frequency bands, we find the 
error to vary between 1 and 3 times the rms error and 
therefore we set it to be 1.5 times the rms error in our 
error estimates. 

Based on all these simulations, we conclude that the 
errors in polarized flux density, EVPA and fractional po- 



larization should be calculated using Eqs. Bl B2 and 



B3 where ctq and au should include contributions from 
the rms error, D-term error and CLEAN error so that 
they are defined as: 



'Dtc 



+ (1.5 X CTrms)^, 



(B5) 



where (Trms is the corresponding rms error and (TDtcrm is 
calculated using Eq. |B4[ The errors for our analysis are 
defined using the above criteria. 

C. ERRORS IN RM MAPS 

Errors in the RM values can formally be obtained from 
the variance-covariance matrix of the linear fit as the er- 
ror of the slope. This method is used in the RM task in 
AIPS. The resulting errors depend largely on the error 
of the EVPA, which is dominated by the absolute cal- 
ibration error in the bright jet locations. As we show 
in Appendix |Bj the rms contribution of the EVPA er- 
ror is well-described by error propagation from Q and U 
rms values. It is crucial to obtain a correct error esti- 
mate for the RM values in order to distinguish between 
various depolarization models and to properly study gra- 
dients in the RM within a source. Using the same kind 
of approach as in Appendix [B] we simulated the effect 
of random noise and calibration error to see if the error 
estimate from the variance-covariance matrix is correct. 

The simulations were done in the same manner as in 
Appendix [B] but repeating the procedure for all the fre- 
quency bands. We then added a random error drawn 
from a Gaussian distribution with standard deviation 
corresponding to the calibration error at the given fre- 
quency band to each of the EVPA images. This error is 
added systematically to each pixel in a given image but 
varies randomly from one epoch to another and needs 
to be accounted for in order to obtain correct error es- 
timates for the RM. The simulated RM maps were then 
created using the same script as for the real data and 
1000 RM maps were obtained. Without any added noise, 
the expected value for the RM in each map is rad m~^ 
minus any Galactic Faraday rotation. This is because we 
created the simulated maps using the exact same proce- 
dure as for the real data and therefore also the Galactic 
Faraday rotation correction was taken into account. 
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Point at 15.3 GHz: rms only 
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Point at 15.3 GHz: D-term noise added 
Expected I = 7.850, <rms> = 0.000450 

Observed <l> = 7.853 c = 0.000845 
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Point at 15.3 GHz: D-term noise added 
Expected Q = 0.555, sigma_0 = 0.00136 

Observed <Q> = 0.56 o = 0.001 03 
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Point at 15.3 GHz: D-term noise added 
Expected U = -0.555, sigma_U = 0.00136 
Observed <U> = -0.56 a = 0.001 25 
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Point at 15.3 GHz: D-term noise added 
Expected EVPA = -22.50, 5igma_EVPA = 0.05 

Observed <EVPA> = -22.5 a = 0.038 
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Point at 15.3 GHz: D-term noise added 
Expected m = 10.0, sigma_m = 0.017 

Observed <m> = 1 0.0 a = 0.01 56 
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Fig. 17. — Distributions of the 100 simulated values of t otal intensity, Stokes Q and U, polarized flux density, EVPA and fractional 
polarization for 1226+023 from jet location shown in Fig. |15| at 15 GHz. Top six panels show simulations including the rms noise only 
and bottom six panels with D-term noise added. The expected value is shown by a dashed line (red in the online version). See text for 
information on the expected error values. 



We repeated the simulations for three sources, 
0333+321, 0735+178 and 1226+023. The source 
0735+178 was included so that we could study the ef- 
fect on a faint source which does not have a well-resolved 
jet. The expected value of RM in each source is differ- 
ent because of the differing Galactic Faraday rotation. 
For 0333+321 we expect a RM of -34.1 rad m-^, for 
0735+178 a RM of -20.5 rad m'^ and for 1226+023 a 
RM of -4.7 rad m^^^ 

C.l. Errors in individual jet locations 

In order to study the RM error in individual jet lo- 
cations, we chose the same jet locations for 0333+321 
and 1226+023 as in Appendix |B] and shown in Figs. [TS] 
and [T6l For 0735+178 we chose the locations shown in 
Fig. |20| The distributions of the 1000 simulated values 
are sEown in Figs. 21-23 Once again it is clear that 



the standard deviations a from our simulations, listed 
below each distribution, agree well with the error from 
the variance-covariance matrix <RMe> (average value 
over the 1000 simulations), listed above each distribu- 
tion. There is, however, a visible offset in the peaks of 
the distributions to the left of the expected value. This 
offset is due to small errors in the final CLEAN proce- 



dure as shown in Appendix B] and it is accounted for by 
the additional error due to the CLEAN procedure, which 
is taken into account in our final Q and U errors. 

Our simulations also show that it is common to obtain 
very high RM values of ^ ±2 x 10^ rad m~^ purely due 
to noise in the data. An example of such simulated jet 
is shown in Fig. [24j We find at least one very high RM 
pixel in 67% of the simulated maps in 0735+178, in 74% 
in 0333+321 and in all the simulated maps of 1226+023. 
The number of the high-RM pixels depends on the total 
number of pixels in the simulated maps. In 1226+023 
the fraction is much higher because it has 5 times more 
pixels than the two other sources. The median fraction 
of high-RM pixels in the simulated maps is 0.3 - 0.6 % 
depending on the source but in 5.5 - 26% of the maps 
there are more than 100 such pixels, resulting in patches 
of at least 1 times 1 mas. This happens because the pix- 
els are not independent on scales smaller than the beam 
size. Therefore even fairly large very high-RM regions 
like these in the real RM maps are most likely spurious. 
For this reason, in our real RM maps, we have blanked 
the high-RM pixels because they are most likely due to 
noise in the data and do not represent real structure. 

Based on the simulations we conclude that the error 



Faraday rotation in parsec-scale AGN jets 



27 



Point 5 at 15.3 GHz: original 
Expected I = 0.0117, <rms> = 0.000407 
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Point 5 at 15.3 GHz: original 
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Point 5 at 15.3 GHz: original 
Expected EVPA = 35.71, <rms> = 10.4 

Observed <EVPA> = 33.7 a = 5.99 
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Point 5 at 15.3 GHz: original 
Expected m = 10.05, <rms> = 3.63 

Observed <m> = 1 4.6 a = 3.25 



Point 5 at 15.3 GHz: clipped and restored 
Expected I = 0.041, <rms> = 0.000434 

Observed <l> = 0.04 c ^ 0.00051 9 
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Point 5 at 15.3 GHz: clipped and restored 
Expected = -0.0013, <rm5> = 0.000455 

Observed <0> = -0.00072 a = 0.000441 
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Observed <P> = 0.0032 a = 0.000447 
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Point 5 at 15.3 GHz: clipped and restored 
Expected EVPA = -54.08, <rms> = 3.2 

Observed <EVPA> = -51 .6 o = 4.05 
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Point 5 at 15.3 GHz: clipped and restored 
Expected m = 10.0, <rms> = 1.13 

Observed <m> = 7.9o-1.14 



Fig. 18. — Distributions of the 100 simulated values of tota l intensity, Stokes Q and U, polarized flux density, EVPA and fractional 
polarization for 0333+321 from jet location 5 shown in Fig. |16| Top six panels show the original 15 GHz data and bottom six panels data 
which has («,f)-range clipped as our real data and which were restored with the 8.1 GHz beam size. The expected value is shown by a 
dashed line (red in the online version). See text for information on the expected error values. 



from the variance-covariance matrix of the hncar fit is a 
good estimate for the true error in the RM in individual 
jet locations. 

C.2. Errors in RM gradients 

Over the past decade there have been several c laims 
of detectio ns of RM gradients in many AGN (e.g 
topoulos et al. 2009, for a compilation of results) 



Con 
hese 

results are still somewhat controversial due to the dif- 



ficulty in addre ss ing the errors in the g r adien ts (Taylor 



k Zavala|2010). Murphy & Gabuzda (2011) convolve 



simulated gradients with ditterent beam sizes and argue 
that a gradient can be significant even when the jet is not 
resolved. Their simulation does not, however, take into 
account noise in the data. With simulations we have 
shown that the errors in individual pixels are well de- 
scribed by the error from the variance-covariance matrix 
of the linear fit. It is more complicated to study the 
errors in the RM gradients, or more precisely, the prob- 
ability of detecting a spurious RM gradient due to noise 
in the data. Several of our simulated jets show struc- 
tures resembling real gradients primarily due to effects 
of the finite beam size. To quantify this effect, we have 
taken several transverse slices across the simulated jets. 



shown in Figs. [25] - [27j and fitted a simple line to the 
data. We have then looked at the distributions of the 
slopes which give us the maximum spurious gradient in 
rad m~^ mas ~^, created by noise in the data. 

In Fig. [28] we show examples of the slices in the dif- 
ferent sources. One thing to note is that when studying 
gradients, the absolute EVPA calibration error can be 
ignored because it affects each pixel in the same direc- 
tion an d thus will not affect a gradient across the same 
source (Mahmud et al. 2009). Therefore in each pixel 
we have subtracted in quadrature the amount of calibra- 
tion error, ~ 60 rad m~^, from the total RM error. It is 
very important to note that the pixels are not indepen- 
dent but very much affected by the beam size. We have 
plotted the size of the beam along the transverse slice in 
each plot to show the scale in which the pixels are not 
independent. As can be seen from the plots, depending 
on the slice, there are only 1-4 independent points along 
each slice. 

T he d istributions of the measured slopes are shown in 
Fig. 
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In the calculation of the slope we have ignored 
the very-high RM values. This resulted in slices with only 
a few pixels with measured RMs which usually resulted 
in steep slopes. These are ignored in the distribution 
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Fig. 19. — Distributions of the 100 simulated values of t ota l intensity, Stokes Q and U, polarized flux density, EVPA and fractional 
polarization for 0333+321 from jet location 5 sliown in Fig. |16| Top six panels show 12.1 GHz data and bottom six panels 8.1 GHz. The 
expected value is shown by a dashed line (red in the online version). See text for information on the expected error values. 



plots. The most notable are the distributions from slices 
C in 0333+321 and slices D in 0735+178 which show gra- 
dients up to 200 rad m~^ mas""'^. The slices are about 
1.5 beam sizes wide and with a typical beam width along 
the slice of 1.5 mas, this would result in a total gradient 
of 450 rad m^^ over the slice. This is of the same order 
as most of the gradients in the literature. It is very close 
to the value we see in the jet of 3C 273. Reassuringly, the 
magnitude of the spurious gradient is strongly dependent 
on the width of the jet with respect to the beam size as 
shown in Fig. 30 left panel. When the width of the po- 
larized jet approaches three beam widths, the spurious 
gradient diminishes. To calculate the jet width in beams 
we first take the width of the slice in pixels and divide 
that by the beam size along the slice direction in pixels. 
This is because we want to ensure that we take into ac- 
count that we do not detect polarization in as large area 
as total intensity and require the jet to be resolved in 
polarization as well. 

We also note that almost all the spurious gradients we 
see in the simulations are smaller than three times the 
RM errors at the location of the slice if using the errors 
as we have defined them and therefore would not be ac- 
cepted as true gradients in our real data. This is shown 



in Fig. [30| right panel, where we plot the fraction of false 
positives, i.e. number of spurious gradients which exceed 
the error limit, against the jet width for la, 2a and 3a 
limits. Here a is defined as the largest error bar in the 
end of a gradient slice. This plot clearly demonstrates 
how the fraction of false positives goes to zero if the jet 
is two beams wide and a 3a error limit is used. When the 
jet is more than 2.5 beams wide, even a 2a limit might 
be sufficient, if additionally taking into account the mag- 
nitude of the maximum possible spurious gradient shown 
in the left panel. For our specific four-frequency setup 
between 8 and 15 GHz using 128 Mbits s~^ recording bit 
rate and ~60 minute on-source time, we find the limit to 
be 200 rad m~^ mas~^. 

Our simulations also verify that it is impossible to get 
a persistent gradient over the whole jet length, as is seen 
in 3C 273, due to noise in the data. Therefore it is im- 
portant that if gradients are detected, they are seen in 
multiple jet locations more than a beamwidth apart, or 
in the same jet location over multiple epochs, and that 
the jet is well-resolved in polarized flux density, prefer- 
ably over two times the beam width. When studying 
jets that are not wide enough, our simulation results can 
be used as "rules of thumb" for estimating the magni- 



Faraday rotation in parsec-scale AGN jets 



29 



0735+178 




-2 



Relative R.A- (mas) 

Fig. 20. — Locations of the simulated RM distributions 
0735+178 overlaid on the 15 GHz total intensity contours. 



tude of a reliable gradient. For example, in a jet that is 
two beams wide, it is possible to have spurious gradients 
up to 100 rad m~^ mas~^ which would result in a total 
gradient of 300 rad m^^ if the beam width is 1.5 mas. 
We note that detection of three 2a gradients at different 
epochs corresponds to a 3cr detection so that multi-epoch 
observations can help to determine if observed gradients 
in jets that are not well-resolved are real. However, this 
requires correct treatment of the RM errors and is not 
recommended for jets less than 1.5 beams wide due to 
the large fraction of false positives at these jet widths. 

Our simulations overall agree well with the suggestions 
by Taylor & Zavala (2010) to consider a gradient reliable 
if the jet is more than three beams wide and use Scr error 
limits. Based on our simulations a jet which is two beams 
wide is already sufhcient if the gradient exceeds the 3cr 
limit. Attempting to detect RM gradients in jets less 
than two beams wide is discouraged as the probability to 
detect false positives exceeds 0.02 even when a 3ct limit 
is used. We emphasize that this approach requires the 
RM errors to be calculated from the variance-covariance 
matrix of the EVPA vs. fit with appropriate EVPA 
errors taken into account; calibration errors can then be 
subtracted from the error bars. These simulations are ap- 
plicable to the four-frequency configuration used in this 
paper and should be repeated if different frequency bands 
are used. 



REFERENCES 



Alberdi, A., Gomez, J.-L., Marcaide, J. M., Marscher, A. P. & 

Perez- Torres, M. A. 2000, A & A, 361, 529 
Algaba, J. C, Gabuzda, D. C. & Smith, P. S. 2011, 

arXiv:1110:5717 
Asada, K., Inoue, M., Uchida, Y. et al. 2002, PASJ, 54, L39 
Asada, K., Inoue, M., Kameno, S., & Nagai, H. 2008, ApJ, 675, 79 
Asada, K., Inoue, M., Nakamura, M., Kameno, S., & Nagai, H. 

2008, ApJ, 682, 798 
Asada, K., Nakamura, M., Inoue, M., Kameno, S., & Nagai, H. 

2010, ApJ, 720, 41 
Attridge, J. M., Wardle, J. F. C. & Homan, D. C. 2005, ApJ, 633, 

85 

Blandford, R,. 1993, in Astrophysical Jets, Astrophysics and 
Space Science Library (Cambridge: Cambridge Univ. Press), 
Vol. 103, 15 

Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433 
Broderick,, A. E. & McKinney, J. C. 2010, ApJ, 725, 750 
Burn, B. J. 1966, MNRAS, 133, 67 
Chen, T. 2005, Ph.D. thesis, Brandeis Univ. 
Clausen-Brown, E., Lyutikov, M. & Kharb, P. 2011, MNRAS, 
415, 2081 

Contopoulos, I., Christodoulou, D. M., Kazanas, D. & Gabuzda, 

D. C. 2009, ApJ, 702, L148 
Croke, S. & Gabuzda, D. C. 2008, MNRAS, 386, 619 
Croke, S. M., O'SuUivan, S. P. & Gabuzda, D. C. 2010, MNRAS, 

402, 259 

Gabuzda, D. C, Murray, E. & Cronin, P. 2004, MNRAS, L89 
Gabuzda, D. C, Vitrishchak, V. M., Mahmud, M. & O'SuUivan, 

S. P. 2008, MNRAS, 384, 1003 
Gomez, J.-L., Marscher, A. P., Alberdi, A., Jorstad, S. G. & 

Agudo, I. 2002, VLBA Scientific Memo No. 30 
Gomez, J.-L., Marscher, A. P., Jorstad, S. G., Agudo, I & 

Roca-Sogorb, M. 2008, ApJ, 681, L69 
Gomez, J.-L., Roca-Sogorb, M., Agudo, I, Marscher, A. P. & 

Jorstad, S. G. 2011, ApJ, 733, 11 
Homan, D. C, Attridge, J. M. & Wardle, J. F. C. 2001, ApJ, 556, 

113 

Homan, D. C. & Lister, M. L. 2006, AJ, 131, 1262 

Homan, D. C, Lister, M. L., Aller, H. D., AUer, M. F. & Wardle, 

J. F. C. 2009, ApJ, 696, 328 
Homan, D. C. 2012, ApJL, 747, L24 



Hovatta, T., Valtaoja, E., Tornikoski, M. & Lahteenmaki, A. 

2009, a, 494, 527 
Hughes, P. A, Aller, H. D. & Aller, M. F. 1989, ApJ, 341, 68 
Jones, T. W. & O'Dell, S. L. 1977, ApJ, 214, 522 
Kellermann, K. I., Lister, M. L., Homan, D. C, et al. 2004, ApJ, 

609, 539 

Komatsu, E., Dunkley J., Nolta, M. R. et al. 2009, ApJS, 180, 
330 

Kovalev, Y.Y., Lobanov, A. P., Pushkarev, A. B. &; Zensus, J. A. 

2008, A&A, 483, 759 
Laing, R. 1980, MNRAS, 193, 439 

Lister, M. L., Aller, H. D., Aller, M. F., et al. 2009a, AJ, 137, 3718 
Lister, M. L., Cohen, M. H., Homan, D. C. et al. 2009b, AJ, 138, 
1874 

Lister, M. L. & Homan, D. C. 2005, AJ, 130, 1389 

Lister, M. L., Tingay, S. J., Murphy, D. W., Finer, B. G., Jones, 

D. L., Preston, R. A. 2001, ApJ, 554, 948 
Lobanov, A. P. 1998, A&A, 330, 79 

Lyutikov, M., Pariev, V. I. & Gabuzda, D. C. 2005, MNRAS, 360, 
869 

Mahmud, M., Gabuzda, D. & Bezrukovs, V. 2009, MNRAS, 400, 2 
Mantovani, F., Rossetti, A., Junor, W., Saikia, D. J. & Salter, C. 

J. 2010, A&A, 518, A33 
Marscher, A. P., Jorstad, S. G., D'Arcangelo, F. D. et al. 2008, 

Nature, 452, 966 
Marr, J. M., Taylor, G. B. & Crawford HI, F. 2001 ApJ, 550, 160 
Meier, D. L., Koide, S. & Uchida, Y. 2001, Science, 291, 84 
McKinney, J. C. & Narayan, R. 2007, MNRAS, 375, 531 
Murphy, E. & Gabuzda, D. C. 2011, arXiv:1109.4778 
Mutel, R. L., Denn, G. R. & Dreier, C. 2005, ASPC, 340, 155 
O'SuUivan, S. P. & Gabuzda, D. C. 2009, MNRAS, 393, 429 
Pacholczyk, A. G. 1970, Radio Astophysics, Nonthermal 

Processes in Galactic and Extragalactic Sources, (W. H. 

Freeman and Company) 
Press, W. H., Flannery, B. P., Teukolsky, S. A., & VetterUng, W. 

T. 1992, Numerical Recipes in FORTRAN 77, 2nd edn. 

(Cambridge University Press) 
Pushkarev, A. B. 2001, ARep, 45, 667 

Pushkarev, A. B., Hovatta, T., Kovalev, Y. Y. et al. 2012, A&A, 

submitted 

Reichstein, A. & Gabuzda, D. 2011, |arXiv:1102.0702| 



30 



Point 0: <l>= 9.64 <RMe>= 70.9 



Hovatta et al. 
Point 1 : <l>= 4.33 <RMe>= 70.9 



Point 2: <l>= 0.68 <RMe>= 71.0 



-400 -200 200 400 

<RM>=-29.5o = 69.7 



-400 -200 200 400 

<RM>=24.6a = 70.2 



-400 -200 200 400 

<RM>=-25.1 = 70.1 



Point 3: <l>= 2.33 <RMe>= 70.9 



Point 4: <l>= 0.28 <RIVIe>= 71 .4 



Point 5: <l>= 0.60 <RMe>= 71.0 




-400 -200 200 400 

<RM>=43.7a = 70.0 



-400 -200 200 400 

<RM>=-63.9o = 71.2 



-400 -200 200 400 

<RIVI>=29.7a = 70.4 



Point 6: <l>= 0.03 <RMe>= 86.5 



Point 7: <l>= 0.98 <RIVIe>= 71.0 



Point 8: <l>= 0.06 <RMe>= 75.7 



J. 



-400 -200 200 400 

<RM>=-85.9o = 78.8 



-400 -200 200 400 

<RIVI>=33.7a = 70.5 



-400 -200 200 400 

<RM>=-33.1 o = 74.2 



Point 9: <l>= 0.05 <RMe>= 77.2 



^111 1^ 



-400 -200 200 400 

<RM>=-27.2a = 76.5 



Fig. 21. — Distributions of the 1000 simulated RM values for 1226+023 from jet locations shown in Fig. |15| Average total intensity value 
<I>, average RM error from the variance-covariance matrix <RMe>, average RM value <RM> from the simulated jets and the standard 
deviation of the distribution are given for each location. The expected RM value is shown by a dashed line (red in the online version). 
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Fig. 22. — Distributions of the 1000 simulated RM values for 0333+321 from jet locations shown in Fig. |16| Average total intensity value 
<I>, average RM error from the variance-covariance matrix <RMe>, average RM value <RM> from the simulated jets and the standard 
deviation of the distribution are given for each location. The expected RM value is shown by a dashed line (red in the online version). 
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Fig. 23. — Distributions of the 1000 simulated RM values for 0735+178 from jet locations shown in Fig. |20| Average total intensity value 
<I>, average RM error from the variance-covariance matrix <RMe>, average RM value <RM> from the simulated jets and the standard 
deviation of the distribution are given for each location. The expected RM value is shown by a dashed line (red in the online version). 
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Fig. 24. — Example of a simulation in which extreme RM values 
are generated due to noise in the data. Color version is available 
in the online edition of the Journal. 
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Fig. 25. — Slices transverse to the jet in 1226+023 overlaid on a 
simulated RM map. Color version is available in the online edition 
of the Journal. 
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Fig. 26. — Slices transverse to the jet in 0333+321 overlaid on a 
simulated RM map. Color version is available in the online edition 
of the Journal. 
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Fig. 27. — Slices transverse to the jet in 0735+178 overlaid on a 
simulated RM map. Color version is available in the online edition 
of the Journal. 
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Fig. 28. — The most extreme examples of simulated transverse gradients over the slices across the jet locations shown in Figs 
None of these simulated jets contain an actual RM gradient. The FWHM beam size along the slice direction is shown in each 
scale bar. 
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Fig. 29. — Distributions of 1000 simulated gradients across the jet locations shown in Figs. [25] - 127| 1226+023 slice A is omitted from 
the plot because it is almost identical to slice B. 




Fig. 30. — Maximum gradient from the simulations against the jet width (left). Fraction of total gradients along the slice in our simulations 
exceeding 1 (filled circles, red in the online version), 2 (filled squares, black in the online version) or 3 (filled diamonds, blue in the online 
version) standard deviations against the jet width. 
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Fig. 31. — Fits of depolarization curves to all isolated component. To appear Online only 
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Fig. 32. — Fits of depolarization curves to all isolated component. To appear Online only 
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Fits of depolarization curves to all isolated component. To appear Online only 
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Fig. 34. — Fits of depolarization curves to all isolated component. To appear Online only 
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Fig. 36. — Fits of depolarization curves to all isolated component. To appear Online only 
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TABLE 4 

Isolated components at 15.3 GHz and their polarization values used in the 

DEPOLARIZATION ANALYSIS 



Source 


I.D. 


Epoch 


r 


P.A. 


RM 


mi5.3 


mi2.i 


in8.4 
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sigiiif. of 








(mas) 


(dog) 


(rad m~^) 


(%) 


(%) 


(%) 


(%) 


(x lO^m"^) 
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(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


0003-066 


4 


2006-07-07 


6.86 


-81.1 


113.5 ± 9.6 


8.6 ± 1.1 


8.2 ± 1.4 


8.2 ± 0.6 


7.9 ± 0.7 


-0.4 ± 0.8 


0.4 


0119+115 


1 


2006-06-15 


0.09 


7.6 


-59.3 ± 8.7 


3.2 ± 0.1 


3.5 ± 0.1 


3.5 ± 0.1 


3.5 ± 0.1 


+0.3 ± 0.1 


2.0"^ 


0133+476 


1 


2006-08-09 


2.89 


-26.6 


-153.5 ± 4.6 


10.8 ± 0.6 


9.7 ± 0.4 


8.4 ± 0.4 


4.4 ± 0.5 


-2.3 ± 0.4 


5.8"= 


0149+218 


1 


2006-02-12 


8.60 


-11.8 




22.7 ± 4.5 


29.1 ± 3.9 


15.1 ± 3.1 


< 10.0 


-4.7 ± 1.7 


2.8 


0212+735 


1 


2006-07-07 


0.02 


-127.4 


-284.9 ± 11.6 


9.0 ± 0.1 


8.5 ± 0.04 


5.4 ± 0.04 


3.6 ± 0.1 


-4.1 ± 0.1 


77.8"= 


0333+321 


7 


2006-07-07 


6.25 


129.7 


-1.2 ± 8.0 


7.0 ± 1.6 


11.7 ± 1.1 


9.6 ± 0.8 


8.7 ± 0.9 


-1.2 ± 0.8 


1.5 


0415+379 


24 


2006-05-24 


7.19 


60.4 




< 11.8 


9.4 ± 2.4 


8.0 ± 1.9 


7.2 ± 1.7 


-2.1 ± 2.1 


1.0 


0429+415 


2 


2006-03-09 


20.54 


59.4 




14.0 ± 1.0 


14.0 ± 0.8 


2.3 ± 0.5 


1.8 ± 0.6 


-11.6 ± 1.5 


7.8"= 


0430+052 


11 


2006-05-24 


9.49 


-117.4 


-2.5 ± 3.4 


7.5 ± 0.5 


8.6 ± 0.4 


7.5 ± 0.3 


7.3 ± 0.3 


-0.7 ± 0.3 


2.1 


0430+052 


13 


2006-05-24 


6.00 


-114.6 


622.0 ± 81.5 


9.0 ± 1.4 


8.0 ± 1.1 


4.3 ± 0.9 


3.9 ± 0.9 


-4.9 ± 1.4 


3.6 


0528+134 


2 


2006-10-06 


3.24 


24.8 


56.8 ± 9.6 


14.7 ± 2.4 


9.7 ± 1.7 


8.4 ± 1.3 


9.9 ± 1.4 


-1.9 ± 1.1 


1.7 


0605-085 


4 


2006-11-10 


2.51 


123.6 


41.4 ± 9.4 


9.3 ± 2.2 


5.5 ± 1.7 


7.4 ± 1.2 


7.4 ± 1.4 


-0.5 ± 1.7 


0.3 


0735+178 


5 


2006-04-28 


3.57 


64.6 




< 5.0 


< 4.0 


4.6 ± 1.2 


5.2 ± 1.3 


+0.7 ± 2.3 


0.3 


0754+100 


2 


2006-04-28 


4.97 


13.8 


-105.1 ± 4.5 


17.7 ± 5.2 


22.3 ± 3.6 


22.5 ± 3.1 


25.4 ± 3.3 


+ 1.0 ± 1.3 


0.8 
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2 


2006-02-12 


3.92 


157.5 


552.1 ± 5.2 


10.1 ± 0.9 


9.5 ± 0.7 


9.2 ± 0.4 


8.0 ± 0.4 


-1.0 ± 0.5 


2.0 
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1.91 
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1.2 ± 0.2 
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-0.7 ± 0.1 


8.4 


1226+023 


12 


2006-06-15 
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3.1 ± 0.1 


2.9 ± 0.04 


0.7 ± 0.04 


0.6 ± 0.1 


-9.6 ± 0.3 


29.5'= 
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4.5 ± 0.5 


3.1 ± 0.9 


5.6 ± 0.3 


6.0 ± 0.3 


+ 1.8 ± 0.8 


2.4 
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5 


2006-04-05 


4.10 


-125.9 




11.2 ± 0.3 


12.7 ± 0.5 


12.0 ± 0.2 


10.7 ± 0.2 


-0.3 ± 0.1 


1.8" 
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1 


2006-09-06 


6.48 


-127.0 
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0.3 
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3 


2006-11-10 


44.67 
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13.3 ± 3.6 
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< 7.4 


-3.5 ± 2.5 


1.4 


1418+546 


7 


2006-02-12 
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4.3 ± 1.3 


-3.0 ± 1.4 


2.1 


1418+546 


7 


2006-11-10 


5.39 


129.3 




< 7.1 


8.0 ± 2.1 


8.9 ± 1.9 


6.8 ± 2.2 


+0.4 ± 2.2 


0.2 


1458+718 


2 


2006-09-06 


27.27 


166.4 




< 19.0 


19.9 ± 3.2 


22.9 ± 2.6 


26.2 ± 3.0 


+ 1.7 ± 1.3 


1.3 


1502+106 


2 


2006-07-07 


2.36 


122.5 




7.1 ± 2.0 


5.3 ± 1.4 


3.3 ± 1.1 
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-3.0 ± 2.3 


1.3 
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1 


2006-03-09 


17.67 
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< 18.9 


< 20.1 


24.4 ± 4.7 


26.4 ± 5.1 
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0.9 
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< 8.6 


8.1 ± 1.9 


21.3 ± 1.6 
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3 
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3.95 
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7.5 ± 1.3 


8.4 ± 0.9 
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7.4 ± 1.3 


+0.9 ± 0.9 


1.0 
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0.3 
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1.4 


1655+077 


2 
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8.17 


-42.1 


25.6 ± 13.5 


19.2 ± 2.0 


19.0 ± 1.8 


13.8 ± 1.3 


14.0 ± 1.41 


-2.1 ± 0.7 


3.0 


1730-130 


2 


2006-07-07 


8.67 


14.1 


-77.2 ± 3.5 


9.5 ± 1.6 


9.3 ± 1.1 


9.3 ± 0.9 


9.9 ± 0.8 


+0.3 ± 0.8 


0.3 


1730-130 


6 


2006-07-07 


3.08 


26.3 


41.6 ± 4.6 


8.9 ± 1.0 


7.4 ± 0.7 


6.4 ± 0.6 


7.4 ± 0.5 


-0.8 ± 0.6 


1.2 


1828+487 


8 


2006-08-09 


2.74 


-38.3 




6.72 ± 1.4 


4.8 ± 1.0 


3.4 ± 0.8 


< 2.70 


-4.2 ± 1.9 


2.2 


1828+487 


5 


2006-08-09 


6.29 


-27.6 


-97.8 ± 15.5 


17.1 ± 2.1 


17.2 ± 1.5 


12.7 ± 1.3 


7.0 ± 1.4 


-3.1 ± 0.9 


3.6 


1828+487 


3 


2006-08-09 


11.58 


-32.5 


-17.3 ± 13.8 


20.6 ± 2.8 


12.7 ± 2.0 


14.4 ± 1.7 


7.7 ± 1.8 


-2.3 ± 1.1 


2.1"= 


1908-201 


1 


2006-03-09 


3.49 


37.6 


-108.8 ± 22.7 


15.8 ± 3.3 


14.4 ± 3.4 


12.2 ± 2.0 


13.6 ± 2.1 


-1.0 ± 1.4 


0.7 


1928+738 


8 


2006-04-28 


2.47 


161.2 


19.0 ± 13.4 


4.5 ± 0.9 


5.9 ± 0.5 


5.9 ± 0.4 


5.6 ± 0.5 


+0.1 ± 0.7 


0.2 


2131-021 


3 


2006-08-09 


1.76 


110.9 


-201.9 ± 15.7 


1.9 ± 0.4 


2.5 ± 0.3 


2.6 ± 0.3 


1.8 ± 0.4 


-0.0 ± 1.0 


0.0 


2155-152 


1 


2006-12-01 


2.54 


-154.7 




4.88 ± 0.4 


3.4 ± 0.3 


2.1 ± 0.5 


3.0 ± 0.4 


-2.7 ± 0.9 


3.0" 


2201+315 


4 


2006-10-06 


3.84 


-140.2 




9.2 ± 2.7 


< 7.6 


8.8 ± 1.9 


7.6 ± 1.9 


-0.4 ± 2.1 


0.2 


2230+114 


2 


2006-02-12 


11.21 


157.6 




< 8.7 


12.4 ± 2.5 


< 5.0 


10.4 ± 1.8 


-1.1 ± 1.7 


0.6 


2251+158 


2 


2006-03-09 


6.08 


-80.1 


-149.8 ± 7.6 


3.5 ± 0.4 


3.6 ± 0.3 


4.3 ± 0.2 


5.1 ± 0.2 


+2.1 ± 0.5 


4.0 


2251+158 
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8.57 


-64.9 


2.8 ± 4.6 


11.7 ± 0.9 


12.8 ± 0.6 


14.5 ± 0.4 


15.2 ± 0.4 


+ 1.2 ± 0.3 


4.1 


2251+158 
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6.07 


-79.9 


26.8 ± 37.5 


2.7 ± 0.4 


4.2 ± 0.3 


4.1 ± 0.2 


3.3 ± 0.3 


+0.1 ± 0.6 


0.2" 


2251+158 


1 


2006-06-15 


8.57 


-64.0 


-124.9 ± 8.0 


16.4 ± 1.0 


14.5 ± 0.7 


16.4 ± 0.5 


12.3 ± 0.7 


-0.1 ± 0.3 


0.4" 



Note. — Columns are as follows: (1) lAU Name (B1950.0); (2) I.D. of the component (0 — core); (3) Epoch; (4) Component distance from the 
phase center of the I map; (5) Position angle of the component from the phase center; (6) RMobs of the component (7) Fractional polarization at 
15.3 GHz. Upper limits arc 3ct limits in fractional polarization; (8) Fractional polarizati on a t 12.1 GHz; (9) Fractional polarization at 8.4 GHz; 
(10) Fractional polarization at 8.1 GHz; (11) Slope b from the depolarization fit in Sect. |4.2| (12) Significance of the slope. Sources marked with 
c have a > 6 indicating low probability that our simple model is accurate in explaining how the polarization changes with wavelength. 
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TABLE 1 

Sources and their rotation measure properties 
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-186.3 


164.5 


1637+574 


OS 562 


0.751 


Q 


12.5 


2006-05-24 


4.3 


294.1 


339.6 


196.3 


1637+826 


NGC 6251 


0.024 


G 




2006-03-09 


-24.3 








1638+398 


NRAO 512 


1.666 


Q 




2006-08-09 


19.4 


-490.0 


-490.0 




1641+399 


3C 345 


0.593 


Q 


19.3 


2006-06-15 


19.4 


59.3 


140.2 


37.0 


1642+690 


4C +69.21 


0.751 


Q 




2006-03-09 


11.5 


81.9 


149.0 


-27.1 


1652+398 


Mrk 501 


0.0337 


B 




2006-02-12 


19.9 


349.3 


395.1 


-46.9 


1655+077 


PKS 1655+077 


0.621 


Q 




2006-11-10 


24.9 


5.6 


-1333.4 


39.1 


1725+044 


PKS 1725+044 


0.293 


Q 


11.7 


2006-03-09 


30.6 


-3846.3 


-3846.3 




1726+455 


S4 1726+45 


0.717 


Q 




2006-09-06 


23.9 
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TABLE 1 — Continued 



Source 


Alias 


z 


Opt. CL 


Papp 

(c) 


Epoch 


Gal. RM 

(rad m^^) 


mcd. RM 
(rad m^^) 


mcd. core RM 
(rad ) 


med. jet RM 
(rad in 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1730-130 


NRAO 530 


0.902 


Q 


27.4 


2006-07-07 


35.9 


15.1 


371.0 


0.9 


1739+522 


4C +51.37 


1.379 


Q 




2006-08-09 


16.3 


-6.1 


-7.0 


28.7 


1741-038 


PKS 1741-03 


1.054 


Q 




2006-12-01 


37.8 


202.7 


204.6 


193.1 


1749+096 


4C +09.57 


0.322 


B 


7.9 


2006-06-15 


64.8 


116.6 


136.6 


63.7 


1749+701 


S4 1749+70 


0.77 


B 




2006-04-05 


8.9 


186.3 


186.3 




1751+288 


B2 1751+28 


1.118 


Q 




2006-10-06 


78.5 


-204.3 


-204.3 




1758+388 




2.092 


Q 




2006-11-10 


53.9 








1800+440 


S4 1800+44 


0.663 


Q 




2006-04-28 


32.7 


244.5 


247.8 


102.0 


1803+784 


S5 1803+784 


0.6797 


B 


10.8 


2006-09-06 


-21.6 


-32.8 


30.0 


-70.8 


1807+698 


3C 371 


0.051 


B 




2006-02-12 


10.8 








1823+568 


4C +56.27 


0.664 


B 


26.2 


2006-07-07 


26.5 


-90.9 


-114.6 


-32.5 


1828+487 


3C 380 


0.692 


Q 


13.1 


2006-08-09 


19.2 


-81.6 


-1172.1 


-66.5 


1845+797 


3C 390.3 


0.0555 


G 




2006-02-12 


-5.3 








1849+670 


S4 1849+67 


0.657 


Q 




2006-05-24 


13.2 


496.7 


496.7 




1901+319 


3C 395 


0.635 


Q 




2006-02-12 


63.7 


-486.9 


-518.1 


107.9 


1908-201 


PKS B1908-201 


1.119 


Q 




2006-03-09 


-20.4 


-202.3 


-365.2 


-130.5 


1928+738 


4C +73.18 


0.302 


Q 


8.2 


2006-04-28 


-20.9 


-863.4 


-962.0 


-3.7 


1936-155 




1.657 


Q 




2006-07-07 


-53.8 


48.1 


48.1 




1957+405 


Cygnus A 


0.0561 


G 




2006-06-15 


-48.2 








1958-179 


PKS 1958-179 


0.652 


Q 




2006-10-06 


-49.5 


-9.8 


-9.8 




2005+403 


TXS 2005+403 


1.736 


Q 


9.9 


2006-09-06 


-40.0 


-141.8 


-642.1 


-108.4 


2008-159 


PKS 2008-159 


1.18 


Q 




2006-11-10 


-38.0 


6456.7 


6495.7 


6401.0 


2021+317 


4C +31.56 




u 




2006-08-09 


-173.0 


-50.8 


-44.7 


-66.5 


2021+614 


OW 637 


0.227 


G 




2006-12-01 


29.1 








2022-077 


PKS 2023-07 


1.388 


Q 




2006-04-05 


-34.9 


-14.5 


-13.5 


-102.7 


2037+511 


3C 418 


1.686 


Q 




2006-05-24 


-4.3 


-467.6 


-502.6 


-197.0 


2043+749 


4C +74.26 


0.104 


Q 




2006-04-28 


-6.5 








2113+293 


B2 2113+29 


1.514 


Q 




2006-02-12 


-95.8 








2121+053 


PKS 2121+053 


1.941 


Q 


11.7 


2006-06-15 


-1.4 


-193.4 


-193.4 




2126-158 




3.28 


Q 




2006-04-28 


0.5 








2128-123 


PKS 2128-12 


0.501 


Q 




2006-10-06 


2.4 


-98.6 


-1316.4 


-95.2 


2131-021 


4C -02.81 


1.284 


Q 




2006-08-09 


4.2 


-34.0 


-31.4 


-129.9 


2134+004 


PKS 2134+004 


1.932 


Q 


5.0 


2006-07-07 


8.9 


374.9 


391.9 


265.4 


2136+141 


OX 161 


2.427 


Q 


4.1 


2006-09-06 


-33.1 


67.4 


67.4 




2141 + 175 


OX 169 


0.2107 


Q 




2006-03-09 


-37.6 








2145+067 


4C +06.69 


0.999 


Q 




2006-10-06 


-14.5 


-197.4 


-350.3 


160.3 


2155-152 


PKS 2155-152 


0.672 


Q 




2006-12-01 


14.2 


-2.6 


-21.7 


23.1 


2200+420 


BL Lac 


0.0686 


B 


10.6 


2006-04-05 


-156.1 


57.0 


-95.8 


62.7 


2200+420 


BL Lac 


0.0686 


B 


10.6 


2006-11-10 


-156.1 


3.7 


41.7 


3.7 


2201 + 171 


PKS 2201+171 


1.076 


Q 




2006-05-24 


-32.1 


114.6 


115.9 


57.3 


2201+315 


4C +31.63 


0.2947 


Q 




2006-10-06 


-120.2 


129.7 


-963.8 


462.4 


2209+236 


PKS 2209+236 


1.125 


Q 




2006-12-01 


-80.5 


68.5 


72.1 


41.0 


2216-038 


PKS 2216-03 


0.901 


Q 




2006-08-09 


0.6 


-22.0 


-130.5 


154.4 


2223-052 


3C 446 


1.404 


Q 


20.3 


2006-10-06 


0.2 


241.7 


241.7 




2227-088 


PHL 5225 


1.5595 


Q 


2.0 


2006-07-07 


-10.9 


-196.7 


-197.8 


-188.3 


2230+114 


CTA 102 


1.037 


Q 


8.6 


2006-02-12 


-25.8 


-4.5 


-190.2 


11.1 


2243-123 


PKS 2243-123 


0.632 


Q 




2006-09-06 


-18.3 


179.8 


456.3 


155.0 


2251 + 158 


3C 454.3 


0.859 


Q 


13.8 


2006-03-09 


-33.5 


-99.2 


-202.3 


-59.4 


2251 + 158 


3C 454.3 


0.859 


Q 


13.8 


2006-06-15 


-33.5 


-77.8 


-64.7 


-80.4 


2254+074 




0.19 


B 


0.6 


2006-04-28 


-11.6 


-4.4 


-7.6 


0.2 


2320-035 


PKS 2320-035 


1.411 


Q 




2006-04-05 


-16.2 


207.8 


395.4 


164.6 


2331+073 


TXS 2331+073 


0.401 


Q 




2006-06-15 


1.4 


591.1 


591.1 




2345-167 


PKS 2345-16 


0.576 


Q 




2006-11-10 


-0.9 


-17.9 


-12.5 


-97.5 


2351+456 


4C +45.51 


1.986 


Q 




2006-05-24 


-27.5 


213.2 


226.9 


48.5 


2356+196 




1.07 


Q 




2006-04-05 


-32.3 








2356+196 




1.07 


Q 




2006-10-06 


-32.3 









Note. — Columns are as follows: (1) lAU Name (B1950.0); (2) alternate name; (3) rcdshift; (4) optical classification where Q — quasar, B — 
BL Lac object, G — active galaxy, and U — unide ntified; (5) Apparen t speed used in viewing angle calculation of Fig. [4| (6) epoch of the RM 
observation; (7) Galactic RM correction taken from [Taylor et al.||2009[ l; (8) median RM over the source; (9) median RM over the core area; (10) 
median RM over the jet. Table 1 is published in its entirety in the electronic edition of the Journal, a portion is shown here for guidance regarding 
its form and content. 
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TABLE 3 



MODELFIT COMPONENTS AT 15.3 CtHz AND THEIR RM VALUES 


Source 


I.D. 


Epoch 


r 


P.A. 


RM 


Isolated 








(mas) 


(deg) 


(rad m~^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0003—066 





2006-07-07 


0.71 


— 168.6 


— 78 ± 73 




OUOo— Odd 


1 


20UO-07-07 


0.66 


— 71.6 


— 88 ± 75 




0003—066 


4 


2006-07-07 


6.86 


—81.1 


113 ± 98 


Y 


0003—066 


5 


2006-07-07 


0.09 


2.9 


—42 ± 72 




0003—066 


6 


2006-07-07 


1.28 


— 102.7 


— 14 ± 77 




0007+106 


1 


2006-06-15 


0.39 


—66.6 


627 ± 130 




0016+731 





2006-08-09 


0.01 


— 52.2 


278 ± 74 




0016+731 


2 


2006-08-09 


0.20 


129.2 


266 ± 74 




0048—097 





2006-09-06 


0.01 


177.8 


— 158 zb 74 




0048—097 


1 


2006-09-06 


1.93 


— 2.3 


— 146 ± 109 




0048—097 


2 


2006-09-06 


0.52 


—8.1 


— 142 ± 75 




0059+581 





2006-05-24 


0.04 


— 17.8 


— 148 ± 72 




0059+581 


3 


2006-05-24 


0.92 


— 115.8 


— 200 ± 77 




0059+581 


4 


2006-05-24 


0.65 


— 154.3 


— 169 ± 73 




0059+581 


5 


2006-05-24 


0.08 


160.7 


— 152 ± 72 




0106+013 





2006-07-07 


0.04 


21.7 


186 ± 88 




0106+013 


2 


2006-07-07 


1.98 


— 105.0 


± 127 




0106+013 


3 


2006-07-07 


1.06 


— 124.4 


— 126 ± 111 




0106+013 


4 


2006-07-07 


0.42 


— 151.7 


151 ± 83 




0119+115 





2006-06-15 


0.16 


— 170.4 


— 88 ± 72 




0119+115 


1 


2()()6-()6-15 


0.09 


7.6 


— 59 ± 72 


Y 


0119+115 


2 


2006-06-15 


1.26 


10.6 


36 ± 74 




0133+476 





2006-08-09 


0.01 


142.7 


— 230 ± 75 




0133+476 


1 


2006-08-09 


2.89 


— 26.6 


— 153 ± 82 


Y 


0133+476 


2 


2006-08-09 


2.72 


— 38.0 


— 129 ± 89 




0133+476 


5 


2006-08-09 


0.48 


—62.5 


— 348 ± 83 




0149+218 


7 


2006-02-12 


0.24 


6.5 


— 1342 ± 74 




0202+319 





2006-08-09 


0.01 


141.0 


— 280 ± 94 




0202+319 


3 


2006-08-09 


0.22 


— 25.6 


— 278 ± 95 




0212+735 





2006-07-07 


1.16 


— 70.9 


—425 ± 79 




0212+735 


1 


2006-07-07 


0.02 


— 127.4 


— 285 ± 71 


Y 


0212+735 


2 


2006-07-07 


0.75 


— 76.8 


— 387 ± 72 




0215+015 





2006-04-28 


0.00 


95.5 


— 102 ± 112 




0215+015 


4 


2006-12-01 


0.23 


125.6 


388 ± 138 




0219+428 


4 


2006-04-05 


2.56 


— 175.3 


180 ± 147 




0219+428 


5 


2006-04-05 


2.08 


— 176.4 


138 ± 111 




0219+428 


5 


2006-11-10 


2.29 


— 177.1 


115 ± 97 




0219+428 


9 


2006-11-10 


1.61 


— 172.8 


48 ± 76 




0219+428 


11 


2006-11-10 


1.00 


— 173.7 


44 ± 72 




0219+428 


13 


2006-11-10 


0.47 


— 172.7 


53 ± 72 




0224+671 





2006-10-06 


0.04 


— 165.6 


— 202 ± 108 




0224+671 


3 


2006-10-06 


0.82 


7.9 


— 68 ± 92 




0224+671 


4 


2006-10-06 


0.49 


10.4 


— 99 ± 90 




0224+671 


8 


2006-10-06 


2.11 


12.6 


54 ± 123 




0234+285 





2006-09-06 


0.36 


163.7 


489 ± 77 




0234+285 


6 


2006-09-06 


0.55 


— 17.1 


202 ± 73 




0234+285 


8 


2006-09-06 


0.21 


8.0 


288 ± 73 




0234+285 


9 


2006-09-06 


0.15 


— 164.0 


388 ± 75 




0235+164 





2006-06-15 


0.00 


108.7 


— 153 ± 80 




0235+164 


1 


2006-06-15 


0.47 


— 15.1 


— 161 ± 80 




0333+321 


7 


2006-07-07 


6.25 


129.7 


— 1 ± 106 


Y 


0333+321 


10 


2006-07-07 


0.99 


128.7 


136 lb 98 




0333+321 


11 


2006-07-07 


1.99 


126.6 


—61 ± 116 




0336—019 





2006-08-09 


0.01 


— 72.3 


151 ± 73 




0336—019 


3 


2006-08-09 


0.39 


77.1 


140 ± 75 




0336—019 


5 


2006-08-09 


1.90 


72.2 


47 ± 73 




0336—019 


6 


2006-08-09 


1.54 


56.9 


99 ± 84 




0336—019 


14 


2006-08-09 


2.21 


53.2 


51 ± 75 




0403—132 





2006-05-24 


0.02 


10.8 


117 ± 78 




0403—132 


4 


2006-05-24 


2.47 


— 171.0 


99 ± 132 




0403—132 


5 


2006-05-24 


0.25 


— 169.0 


122 ± 77 




0414—189 





2006-04-05 


0.00 


— 127.8 


142 ± 146 




0422+004 





2006-06-15 


0.01 


—31.8 


— 13 ± 72 




0422+004 


2 


2006-06-15 


0.71 


4.4 


—3 ± 73 




0429+415 


5 


2006-03-09 


1.00 


9.1 


1771 ± 71 




0429+415 


6 


2006-03-09 


0.09 


— 155.5 


1977 ± 71 




n/l 9Q_L/1 1 ^ 


7 


zuuD-uo-uy 


1 .U 1 


1 fin n 


ori/l S -t- 79 




0430+052 


11 


2006-05-24 


9.49 


-117.4 


-2 ± 77 


Y 


0430+052 


12 


2006-05-24 


8.82 


-118.6 


63 ± 82 




0430+052 


13 


2006-05-24 


6.00 


-114.6 


622 ± 132 


Y 


0430+289 





2006-04-28 


0.00 


127.3 


29 ± 96 




0430+289 


1 


2006-04-28 


0.72 


47.1 


50 ± 125 




0440-003 





2006-07-07 


0.06 


39.9 


110 ± 79 




0440-003 


1 


2006-07-07 


1.26 


-126.8 


-97 ± 86 




0440-003 


2 


2006-07-07 


0.57 


-142.1 


84 ± 75 




0446+112 





2006-09-06 


0.00 


-48.5 


-516 ± 113 




0446+112 


2 


2006-09-06 


0.86 


96.8 


178 ± 124 





Hovatta et al. 

TABLE 3 — Continued 



Source 


I.D. 


Epoch 


r 


P.A. 


RM 












fnias) 


(dog) 


(rad m"^) 


(1) 


[-) 


( 


■> 

■V 




U) 


{'>) 




0458-020 





2006- 


■11- 


■10 


0.00 


-24.0 


-315 ± 80 


0528+134 





2006- 


■10- 


■06 


0.36 


-107.7 


130 ± 71 


0528+134 


2 


2006- 


■10- 


■06 


3.24 


24.8 


57 ± 116 


0528+134 


3 


2006- 


■10- 


■06 


1.77 


36.4 


158 ± 94 


0528+134 


6 


2006- 


-10- 


■06 


0.21 


28.7 


158 ± 71 


0528+134 


7 


2006- 


■10- 


■06 


0.07 


137.5 


150 ± 71 


0529+075 





2006- 


■08- 


■09 


0.07 


-59.2 


128 ± 76 


0529+075 


2 


2006- 


■08- 


■09 


0.66 


-27.2 


138 ± 77 


0529+483 





2006- 


■05- 


■24 


0.08 


-175.2 


183 ± 91 


0529+483 


5 


2006- 


■05- 


■24 


0.17 


6.2 


206 ± 96 


0552+398 





2006- 


■07- 


■07 


0.00 


157.2 


246 ± 75 


0552+398 


1 


2006- 


■07- 


■07 


0.63 


-74.2 


203 ± 91 


0605-085 


4 


2006- 


■11- 


■10 


2.51 


123.6 


41 ± 148 


0605-085 


6 


2006- 


■11- 


■10 


0.72 


130.2 


-138 ± 87 


0607—157 





2006- 


■09- 


■06 


0.01 


-60.0 


— 221 ± 76 


0607-157 


3 


2006- 


■09- 


■06 


0.49 


63.1 


-241 ± 73 


0642+449 


3 


2006- 


■10- 


■06 


0.54 


88.6 


-464 ± 79 


0707+476 





2006- 


■04- 


■05 


0.03 


157.3 


142 ± 96 


0707+476 


7 


2006- 


■04- 


■05 


0.86 


-16.1 


88 ± 106 


0707+476 


8 


2006- 


■04- 


■05 


0.16 


-20.8 


103 ± 94 


0716+714 





2006- 


■0.5- 


■24 


0.01 


-159.6 


-134 ± 81 


0716+714 


6 


2006- 


■05- 


■24 


0.36 


22.5 


-142 ± 82 


0716+714 





2006- 


■12- 


■01 


0.02 


-20.5 


-70 ± 81 


0716+714 


5 


2006- 


■12- 


■01 


1.13 


21.6 


-98 ± 116 


0716+714 


6 


2006- 


■12- 


■01 


0.66 


10.1 


-95 ± 83 


0716+714 


7 


2006- 


■12- 


■01 


0.36 


3.1 


-73 ± 81 


0723-008 





2006- 


■07- 


■07 


0.37 


130.7 


1434 ± 86 


0723-008 


5 


2006- 


■07- 


■07 


0.50 


-48.7 


1515 ± 81 


0723-008 


6 


2006- 


■07- 


■07 


0.20 


-119.5 


1525 ± 78 


0723-008 


7 


2006- 


■07- 


■07 


0.07 


9.8 


1498 ± 79 


0727-115 





2006- 


■10- 


■06 


0.00 


-15.2 


23 ± 72 


0727-115 


3 


2006- 


■10- 


■06 


0.36 


-91.5 


102 ± 72 


0727-115 


4 


2006- 


■10- 


■06 


0.65 


-77.8 


132 ± 72 


0735+178 





2006- 


■04- 


■28 


0.20 


98.3 


35 ± 85 


0735+178 


1 


2006- 


■04- 


■28 


0.66 


73.0 


56 ± 90 


0735+178 


2 


2006- 


■04- 


■28 


1.26 


66.2 


-68 ± 133 


0736+017 





2006- 


■06- 


■15 


0.00 


113.1 


-617 ± 106 


0738+313 


3 


2006- 


■09- 


■06 


4.22 


152.0 


-7 ± 107 


0738+313 


4 


2006- 


■09- 


■06 


3.85 


167.0 


-232 ± 95 


0738+313 


6 


2006- 


■09- 


■06 


3.53 


174.7 


-358 ± 84 


0738+313 


7 


2006- 


■09- 


■06 


2.90 


176.2 


-495 ± 83 


0738+313 


8 


2006- 


■09- 


■06 


1.67 


-174.3 


-728 ± 80 


0748+126 





2006- 


■08- 


■09 


0.01 


-97.7 


-264 ± 73 


0748+126 


4 


2006- 


■08- 


■09 


0.84 


77.2 


-248 ± 80 


0748+126 


6 


2006- 


■08- 


■09 


0.22 


90.5 


-275 ± 73 


0754+100 





2006- 


■04- 


■28 


0.02 


— 141.5 


8 ± 76 


0754+100 


2 


2006- 


■04- 


■28 


4.97 


13.8 


-105 ± 134 


0754+100 


4 


2006- 


■04- 


■28 


1.35 


20.2 


-4 ± 72 


0754+100 


5 


2006- 


■04- 


■28 


0.77 


17.5 


-55 ± 74 


0754+100 


6 


2006- 


■04- 


■28 


1.09 


34.9 


16 ± 72 


0804+499 





2006- 


■10- 


■06 


0.01 


93.9 


-63 ± 106 


0808+019 





2006- 


■08- 


■09 


0.05 


-22.7 


277 ± 72 


0808+019 


1 


2006- 


■08- 


■09 


1.82 


-170.1 


353 ± 137 


0808+019 


2 


2006- 


■08- 


■09 


0.84 


179.8 


295 ± 74 


0808+019 


5 


2006- 


■08- 


■09 


0.16 


155.7 


274 ± 72 


0814+425 





2006- 


■11- 


■10 


0.01 


-95.6 


172 ± 100 


0814+425 


2 


2006- 


■11- 


■10 


1.58 


89.3 


48 ± 119 


0814+425 


3 


2006- 


■11- 


■10 


1.14 


109.4 


—42 ± 101 


0814+425 


6 


2006- 


■11- 


■10 


0.09 


74.0 


176 ± 100 


0823+033 





2006- 


■06- 


■15 


0.00 


-145.6 


-123 ± 73 


0823+033 


8 


2006- 


■06- 


■15 


1.14 


31.2 


-124 ± 87 


0823+033 


9 


2006- 


■06- 


■15 


0.51 


35.9 


-120 ± 73 


0827+243 





2006- 


■05- 


■24 


0.00 


91.0 


137 ± 74 


0827+243 


4 


2006- 


■05- 


■24 


0.87 


100.8 


31 ± 90 


0829+046 





2006- 


■07- 


■07 


0.01 


-53.3 


-147 ± 82 


0829+046 


4 


2006- 


■07- 


■07 


0.61 


62.1 


-12 ± 76 


0829+046 


7 


2006- 


■07- 


■07 


1.44 


63.3 


57 ± 111 


0834-201 





2006- 


■03- 


■09 


0.08 


91.2 


230 ± 77 


Ooo4— 201 


2 


2006- 


■03- 


■09 


0.68 


— 149.7 


35 ± 84 


0834-201 


3 


2006- 


■03- 


■09 


0.04 


-67.2 


230 ± 78 


0834-201 


2 


2006- 


■12- 


■01 


0.62 


-140.1 


102 ± 82 


0834-201 


3 


2006- 


■12- 


■01 


0.06 


-65.1 


143 ± 100 


0836+710 


4 


2006- 


■09- 


■06 


2.94 


-144.6 


-97 ± 114 


0847-120 





2006- 


■02- 


■12 


0.01 


-84.1 


512 ± 90 


0847-120 


3 


2006- 


■02- 


■12 


0.25 


85.9 


483 ± 98 


0847-120 





2006- 


■12- 


■01 


0.05 


163.6 


393 ± 134 


0851+202 





2006- 


■04- 


■28 


0.01 


-167.4 


-367 ± 71 


0851+202 


9 


2006- 


■04- 


■28 


1.58 


-117.6 


20 ± 106 



Isolated 



Faraday rotation in parsec-scale AGN jets 



TABLE 3 — Continued 



Source 


I.D. 


Epoch 


r 
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P.A. 

(dog) 
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(rad m~^) 


Isolated 
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lU 
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Q 
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04 
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2 
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Y 
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u 
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49 


1 q^ ^ 
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1 09 -1- 78 
IZZ Hz 1 o 
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5 
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1 q 8 -1- 77 
loo zL 1 1 
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D 


9nnfi 1 n nR 

ZUUD-IU-UD 


n An 

U.4U 
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uy 1 ( -\-o^^ 


u 
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u.uy 
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lOD.O 
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— ooD ziz yz 




HQ! 7 1 (iOA 

uyj- 1 -\-Oji^ 


D 


onnR na no 
zuuD-uo-uy 


n y18 


ACk ^ 

— 4y .o 


981 -1- 1 n9 
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uyj- / 


/ 


onnft na no 
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12 


K'i R 
— OO.D 


qi n -1- on 
— oiu ziz yu 




uyzo-i-oyz 
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onnfi n7 n7 

ZUUD-U i-V 1 


45 


aq R 

OO.D 
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— IDO i ( 1 




uyzo-i-oyz 
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ZUUD-U i-yj 1 


72 
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— y u i ( z 
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9nnfi n7 n7 
ZUUD-U i-yj 1 


08 
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9 
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1.10 
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— 42 ± 71 




1253-055 


10 


2006-04-05 


0.73 


-142.3 


-70 ± 71 




1253-055 


1 


2006-09-06 


6.48 


-127.0 


397 ± 111 


Y 


1253-055 


4 


2006-09-06 


0.30 


-124.1 


-303 ± 71 




1253-055 


10 


2006-09-06 


0.93 


-140.7 


-101 ± 71 




1308+326 


8 


2006-07-07 


0.76 


-39.2 


-161 ± 77 




1308+326 


11 


2006-07-07 


3.05 


-41.4 


-206 ± 110 




1308+326 


13 


2006-07-07 


3.36 


-46.7 


-254 ± 101 




1308+326 


15 


2006-07-07 


0.31 


-40.5 


-167 ± 81 




1406-076 


7 


2006-04-05 


0.62 


-88.9 


307 ± 110 
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Source I.D. Epoch r P.A. RM Isolated 

(mas) fdcp;) (rad ni^") 



1 1 \ 
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( 






1 1 1 
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-64.8 


415 ± 83 
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2 


2006- 
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4 
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-65 ± 98 




1508-055 
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1 ± 76 
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0.03 
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1510-089 


8 
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-33.2 


63 ± 91 
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9 


2006- 
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-31.7 
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1538+149 
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2 
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3 
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4 
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Y 


1546+027 
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■08- 
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0.16 


-5.9 


53 ± 72 




1546+027 


6 


2006- 


■08- 


■09 


0.48 


172.5 


64 ± 72 




1548+056 





2006- 


■ 10- 


■06 


1.09 


166.3 


-21 ± 72 




1548+056 


3 


2006- 


■10- 


■06 


1.52 


22.5 


-90 ± 82 




1548+056 


4 


2006- 


■10- 


■06 


0.37 


64.5 


-33 ± 71 




1548+056 


5 


2006- 


■10- 


■06 


0.96 


-12.0 


-13 ± 71 




1548+056 


8 


2006- 


■10- 


■06 


0.02 


-169.5 


-28 ± 71 




1606+106 





2006- 


■07- 


■07 


0.00 


140.0 


322 ± 105 




1606+106 


8 


2006- 


■07- 


■07 


1.34 


-57.2 


155 ± 116 




1606+106 


9 


2006- 


■07- 


■07 


0.52 


-44.3 


256 ± 107 




1611+343 


3 


2006- 


■06- 


■15 


3.95 


-168.5 


54 ± 110 


Y 


1611+343 


4 


2006- 


■06- 


■15 


3.07 
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-27 ± 74 




1611+343 


10 


2006- 
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1.64 
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15 ± 78 
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13 
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4.61 


170.4 
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10 
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18 


2006- 


■09- 


■06 


3.38 
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9 
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382 ± 107 
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10 
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0.59 
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391 ± 118 
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0.04 


131.7 
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2 
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■11- 
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8.17 


-42.1 


26 ± 94 


Y 


1655+077 


9 


2006- 


■11- 


■10 


0.48 


-40.1 


-1465 ± 93 




1655+077 


10 


2006- 


■11- 


■10 


0.21 


-44.4 


-1355 ± 90 




1725+044 





2006- 


■03- 


■09 


0.01 


-86.4 


-3850 ± 118 




1730-130 





2006- 


■07- 


■07 


0.04 


-174.4 


515 ± 95 




1730-130 


2 


2006- 


■07- 


■07 


8.67 


14.1 


-77 ± 102 


Y 


1730-130 


6 


2006- 


■07- 


■07 


3.08 


26.3 


42 ± 91 


Y 


1730-130 


8 


2006- 


■07- 


■07 


0.68 


-3.8 


493 ± 94 




1739+522 





2006- 


■08- 


■09 


0.07 


-25.3 


-25 ± 74 




1739+522 


3 


2006- 


■08- 


■09 


0.15 


143.6 


-24 ± 75 




1741-038 





2006- 


■12- 


■01 


0.21 


-174.2 


218 ± 72 




1741-038 


1 


2006- 


■12- 


■01 


0.69 


-142.2 


270 ± 73 




1749+096 





2006- 


■06- 


■15 


0.00 


-109.1 


114 ± 73 




1749+096 


9 


2006- 


■06- 


■15 


0.59 


15.3 


77 ± 72 




1751+288 





2006- 


■10- 


■06 


0.03 


168.2 


-204 ± 97 




1751+288 


1 


2006- 


■10- 


■06 


0.36 


3.4 


-239 ± 102 





Faraday rotation in parsec-scale AGN jets 

TABLE 3 — Continued 



Source 


I.D. 


Epoch 


r 


P.A. 


RM 








fnias) 


(dog) 


(rad m"^) 


(1) 


[-) 




U) 
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1800+440 





2006-04-28 


0.01 


16.0 


289 ± 83 


1800+440 


10 


2006-04-28 


1.09 


-148.9 


105 ± 134 


1800+440 


11 


2006-04-28 


0.24 


-165.9 


265 ± 82 


1803+784 





2006-09-06 


0.01 


104.3 


112 ± 74 


1803+784 


3 


2006-09-06 


1.63 


-93.4 


-51 ± 79 


1803+784 


4 


2006-09-06 


1.38 


-84.1 


-52 ± 76 


1803+784 


5 


2006-09-06 


0.78 


-75.4 


-21 ± 73 


1803+784 


6 


2006-09-06 


0.30 


-81.0 


49 ± 73 


1823+568 





2006-07-07 


0.01 


45.1 


-143 ± 72 


1823+568 


8 


2006-07-07 


1.73 


-161.8 


-17 ± 102 


1823+568 


9 


2006-07-07 


0.40 


-160.5 


-107 ± 72 


1823+568 


10 


2006-07-07 


0.83 


-161.0 


-74 ± 72 


1828+487 


3 


2006-08-09 


11.58 


-32.5 


-17 ± 120 


1828+487 


5 


2006-08-09 


6.29 


-27.6 


-98 ± 107 


1908-201 





2006-03-09 


0.18 


-175.9 


-492 ± 76 


1908-201 


1 


2006-03-09 


3.49 


37.6 


-109 ± 136 


1908-201 


2 


2006-03-09 


0.15 


4.5 


-446 ± 76 


1908-201 


4 


2006-03-09 


2.04 


6.6 


-164 ± 133 


1928+738 


8 


2006-04-28 


2.47 


161.2 


19 ± 99 


1928+738 


9 


2006-04-28 


1.90 


146.5 


39 ± 90 


1928+738 


11 


2006-04-28 


0.12 


150.1 


-1158 ± 73 


1928+738 


12 


2006-04-28 


0.15 


-20.9 


-1051 ± 73 


1928+738 


13 


2006-04-28 


0.48 


-25.3 


-917 ± 75 


1958-179 





2006-10-06 


0.01 


-123.4 


-9 ± 73 


1958-179 


3 


2006-10-06 


0.33 


124.0 


-2 ± 73 


2005+403 





2006-09-06 


0.08 


-59.7 


-753 ± 89 


2005+403 


3 


2006-09-06 


2.12 


88.2 


-81 ± 84 


2005+403 


4 


2006-09-06 


2.12 


99.7 


-130 ± 80 


2005+403 


5 


2006-09-06 


1.57 


102.2 


45 ± 101 


2005+403 


6 


2006-09-06 


0.34 


109.0 


-665 ± 91 


2021+317 





2006-08-09 


0.45 


-37.5 


-31 ± 86 


2021+317 


9 


2006-08-09 


0.44 


176.3 


-69 ± 79 


2021+317 


10 


2006-08-09 


0.14 


135.4 


-65 ± 78 


2021+317 


11 


2006-08-09 


0.18 


-24.1 


-52 ± 80 


2022-077 





2006-04-05 


0.00 


-141.7 


-24 ± 89 


2037+511 


5 


2006-05-24 


0.85 


-147.8 


-543 ± 99 


2037+511 


6 


2006-05-24 


0.35 


-148.2 


-471 ± 78 


2121+053 





2006-06-15 


0.14 


49.7 


-183 ± 72 


2121+053 


3 


2006-06-15 


0.18 


-138.2 


-191 ± 72 


2121+053 


5 


2006-06-15 


0.37 


-105.8 


-199 ± 72 


2128—123 


3 


2006-10-06 


4.67 


-148.6 


-115 ± 110 


2128-123 


4 


2006-10-06 


3.82 


-148.5 


-91 ± 126 


2131-021 





2006-08-09 


0.04 


-80.5 


4 ± 72 


2131-021 


3 


2006-08-09 


1.76 


110.9 


-202 ± 117 


2131-021 


4 


2006-08-09 


0.95 


92.0 


-34 ± 77 


2131-021 


5 


2006-08-09 


0.30 


92.1 


-33 ± 72 


2134+004 


1 


2006-07-07 


1.99 


-105.8 


301 ± 71 


2134+004 


2 


2006-07-07 


1.42 


-89.0 


373 ± 72 


2134+004 


3 


2006-07-07 


1.59 


-119.7 


365 ± 71 


2134+004 


4 


2006-07-07 


0.73 


-109.5 


387 ± 72 


2136+141 


2 


2006-09-06 


1.34 


-130.0 


90 ± 109 


2136+141 


3 


2006-09-06 


0.86 


-105.3 


112 ± 79 


2155-152 





2006-12-01 


0.88 


40.9 


53 ± 82 


2155-152 


3 


2006-12-01 


0.40 


-138.0 


-26 ± 71 


2155-152 


4 


2006-12-01 


0.01 


36.1 


-34 ± 72 


2200+420 
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2006-04-05 


4.23 


-165.4 


53 ± 74 


2200+420 


13 


2006-04-05 


3.22 


-164.2 


66 ± 72 


2200+420 
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2006-04-05 


1.94 


177.8 


48 ± 72 


2200+420 
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1.39 
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2200+420 
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0.34 


82.1 


247 ± 71 


2227-088 





2006-07-07 


0.02 


155.3 


-183 ± 74 
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Note. — Columns are as follows: (1) lAU Name (B1950.0); (2) I.D. of the component (0 = core); (3) Epoch; (4) Component distance from the 
phase center of the I map; (5) Position angle of the component from the phase center; (6) Component RM; (7) Flag for isolated jet components. 
Table 3 is published in its entirety in the electronic edition of the Journal; a portion is shown here for guidance regarding its form and content. 



